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Abstract 

These notes aim to provide a concise pedagogical introduction to some important applications of the renormaliza- 
tion group in statistical physics. After briefly reviewing the scaling approach and Ginzburg-Landau theory for criti- 
cal phenomena near continuous phase transitions in thermal equilibrium, Wilson's momentum shell renormalization 
group method is presented, and the critical exponents for the scalar cD 4 model are determined to first order in a dimen- 
sional e expansion about the upper critical dimension d c = 4. Subsequently, the physically equivalent but technically 
more versatile field-theoretic formulation of the perturbational renormalization group for static critical phenomena is 
described. It is explained how the emergence of scale invariance connects ultraviolet divergences to infrared singular- 
ities, and the renormalization group equation is employed to compute the critical exponents for the (9(n)-symmetric 
Landau-Ginzburg-Wilson theory to lowest non-trivial order in the e expansion. The second part of this overview is 
devoted to field theory representations of non-linear stochastic dynamical systems, and the application of renormaliza- 
tion group tools to critical dynamics. Dynamic critical phenomena in systems near equilibrium are efficiently captured 
through Langevin stochastic equations of motion, and their mapping onto the Janssen-De Dominicis response func- 
tional, as exemplified by the field-theoretic treatment of purely relaxational models with non-conserved (model A) and 
conserved order parameter (model B). As examples for other universality classes, the Langevin description and scaling 
exponents for isotropic ferromagnets at the critical point (model J) and for driven diffusive non-equilibrium systems 
are discussed. Finally, an outlook is presented to scale-invariant phenomena and non-equilibrium phase transitions 
in interacting particle systems. It is shown how the stochastic master equation associated with chemical reactions or 
population dynamics models can be mapped onto imaginary-time, non-Hermitian "quantum" mechanics. In the con- 
tinuum limit, this Doi-Peliti Hamiltonian is in turn represented through a coherent-state path integral action, which 
allows an efficient and powerful renormalization group analysis of, e.g., diffusion-limited annihilation processes, and 
of phase transitions from active to inactive, absorbing states. 
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1. Introduction 

Since Ken Wilson's seminal work in the early 1970s 
1 1 1, based also on the groundbreaking foundations laid 
by Leo Kadanoff, Ben Widom, Michael Fisher (2), 
and others in the preceding decade, the renormalization 



group (RG) has had a profound impact on modern statis- 
tical physics. Not only do renormalization group meth- 
ods provide a powerful tool to analytically describe and 
quantitatively capture both static and dynamic critical 
phenomena near continuous phase transitions that are 
governed by strong interactions, fluctuations, and cor- 
relations, they also allow us to address physical prop- 
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erties associated with the emerging generic scale invari- 
ance in certain entire thermodynamic phases, many non- 
equilibrium steady states, and in relaxation phenomena 
towards either equilibrium or non-equilibrium station- 
ary states. In fact, the renormalization group presents 
us with a conceptual framework and mathematical lan- 
guage that has become ubiquitous in the theoretical de- 
scription of many complex interacting many-particle 
systems encountered in nature. One may even argue 
that the fundamental RG notions of universality and rel- 
evance or irrelevance of interactions and perturbations, 
and the accompanying systematic coarse-graining pro- 
cedures are of crucial importance for any attempt at cap- 
turing natural phenomena in terms of only a few meso- 
or macroscopic degrees of freedom, and thus also form 
the essential philosophical basis for any computational 
modeling, including Monte Carlo simulations. 

In these lecture notes, I aim to give a pedagogical 
introduction and concise overview of first the classic 
applications of renormalization group methods to equi- 
librium critical phenomena, and subsequently to the 
study of critical dynamics, both near and far away from 
thermal equilibrium. The second half of this article 
will specifically explain how the stochastic dynamics of 
interacting many-particle systems, mathematically de- 
scribed either through (coupled) non-linear Langevin or 
more "microscopic" master equations, can be mapped 
onto dynamical field theory representations, and then 
analyzed by means of RG-improved perturbative expan- 
sions. In addition, it will be demonstrated how exploit- 
ing the general structure of the RG flow equations, fixed 
point conditions, and prevalent symmetries yields cer- 
tain exact statements. Other authors contributing to this 
volume will discuss additional applications of renormal- 
ization group tools to a broad variety of physical sys- 
tems and problems, and also cover more recently devel- 
oped efficient non-perturbative approaches. 

2. Critical Phenomena 

We begin with a quick review of Landau's generic 
mean-field treatment of continuous phase transitions in 
thermal equilibrium, define the critical exponents that 
characterize thermodynamic singularities, and then ven- 
ture to an even more general description of critical phe- 
nomena by means of scaling theory. Next we generalize 
to spatially inhomogeneous configurations, investigate 
critical infrared singularities in the two-point correlation 
function, and analyze the Gaussian fluctuations for the 
ensuing Landau-Ginzburg- Wilson Hamiltonian (scalar 
Euclidean O 4 field theory). This allows us to identify 
d c — 4 as the upper critical dimension below which 



2 

fluctuations crucially impact the critical power laws. Fi- 
nally, we introduce Wilson's momentum shell renormal- 
ization group approach, reconsider the Gaussian model, 
discuss the general emerging structure, and at last per- 
turbatively compute the fluctuation corrections to the 
critical exponents to first order in the dimensional ex- 
pansion parameter e = d c -d. Far more detailed exposi- 
tions of the contents of this chapter can be found in the 
excellent textbooks [3|-[8| and in chap. 1 of Ref. fl9). 

2.1. Continuous phase transitions 

Different thermodynamic phases are characterized by 
certain macroscopic, usually extensive state variables 
called order parameters; examples are the magnetization 
in ferromagnetic systems, polarization in ferroelectrics, 
and the macroscopically occupied ground-state wave 
function for superfluids and superconductors. We shall 
henceforth set our order parameter to vanish in the high- 
temperature disordered phase, and to assume a finite 
value in the low-temperature ordered phase. Landau's 
basic construction of a general mean-field description 
for phase transitions relies on an expansion of the free 
energy (density) in terms of the order parameter, natu- 
rally constrained by the symmetries of the physical sys- 
tem under consideration. For example, consider a scalar 
order parameter (f> with discrete inversion or Z2 sym- 
metry that in the ordered phase may take either of two 
degenerate values <p ± = ±\<f>o\. We shall see that the fol- 
lowing generic expansion (with real coefficients) indeed 
describes a continuous or second-order phase transition: 

m= r -4> 2 + ^<t> 4 + ...-hct>, (i) 

if the temperature-dependent parameter r changes sign 
at T c , For simplicity, and again in the spirit of a regular 
Taylor expansion, we let r = a(T - T®), where T® de- 
notes the mean-field critical temperature. Stability re- 
quires that u > (otherwise more expansion terms need 
to be added); near the critical point we can simply take 
u to be a constant. Note that the external field h, thermo- 
dynamically conjugate to the order parameter, explicitly 
breaks the assumed Z2 symmetry (f> — > -(p. 

Minimizing the free energy with respect to <p then 
yields the thermodynamic ground state. Thus, from 
/'((/>) = 0we immediately infer the equation of state 

h(T,cf>) = r(T)<f>+^4>\ (2) 
o 

and the minimization or stability condition reads < 
/"(0) = r + § (f> 2 . At T = T°, {2} reduces to the critical 
isotherm h(T^,(p) — | (fr 1 . For r > 0, the spontaneous 
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order parameter at zero external field h — vanishes; 
for r < 0, one obtains <j> ± = +</>o, where 



b = (6\r\/u) 



1/2 



(3) 



Note the emergence of characteristic power laws in the 
thermodynamic functions that describe the properties 
near the critical point located at T — T®, h = 0. 

The continuous, but non-analytic onset of sponta- 
neous ordering is the hallmark of a second-order phase 
transition, and induces additional thermodynamic sin- 
gularities at the critical point: The isothermal order pa- 
rameter susceptibility becomes Vx T l = (dhjdfi)j - 
r + | 0{L whence 



Xt_ 
V 



l/r 
l/2\r\ 



r>0 
r<0 



diverging as \T - T^\ _1 on both sides of the phase transi- 
tion, with amplitude ratio \t (J I T®)/xt(T T T®) = 2. 
Inserting (OJ into the Landau free energy ([TJ one finds 
for T < Tfand h = 



f(<P ± ) = 



- J* 
4 TU ~ 2u ' 

and consequently for the specific heat 



(d 2 f\ 3a 2 
C h=0 = -VT[-=L\ =VT — 
\dT 2 ) h=Q u 



(5) 



(6) 



whereas per construction /(0) = and C/,=o = in 
the disordered phase. Thus, Landau's mean-field theory 
predicts a critical point discontinuity AC/, = o = VT^ 
for the specific heat. Experimentally, one indeed ob- 
serves singularities in thermodynamic observables and 
power laws at continuous phase transitions, but often 
with critical exponents that differ from the above mean- 
field predictions. Indeed, the divergence of the order 
parameter susceptibility (|4]) indicates violent fluctua- 
tions, inconsistent with any mean-field description that 
entirely neglects such fluctuations and correlations. 

2.2. Scaling theory 

The emergence of scale-free power laws suggests the 
following general scaling hypothesis for the free en- 
ergy, namely that its singular contributions near a criti- 
cal point (T = T c , h = 0) can be written as a generalized 
homogeneous function 



/sing(T,A) 



-\ 2 - a f±{h/\T\ A ) 



(7) 



where r = —f^ measures the deviation from the (true) 
critical temperature T c . Thus, the free energy near criti- 
cality is not an independent function of the two intensive 



control parameters T or t and h, but satisfies a remark- 
able two-parameter scaling law, with analytic scaling 
functions f±(x) respectively for T > T c and T < T c 
that only depend on the ratio x = h/\r\ A , and satisfy 
/±(0) = const. In Landau theory, the corresponding crit- 
ical exponents are a — 0, compare (B), and A = 3/2, as 
can be inferred by combining Eqs. (|2ji and ([3]). The as- 
sociated specific heat singularity follows again via 



Ch=< 



h=0 



VT /3 2 /smg 



T} <9t 2 



C ± \TV a 



(8) 



/!=() 



indicating a divergence if a > 0, and a cusp singularity 
for a < 0. Similarly, one obtains the equation of state 

(4) ^(t, h) = - (^p) = -iri 2 — a /: (h/\ T \ & ) , (9) 



and therefrom the coexistence line at h — 0, t < 



(T,0) 



f-(0) oc \rf 



(10) 



where we have identified /3 = 2 - a - A. 

Additional scaling relations, namely identities that 
relate different critical exponents, can be easily derived; 
for example, on the critical isotherm at t = 0, the in- 
dependence in f' ± on the r.h.s. of |9j must cancel the 
singular prefactor, i.e., f' ± (x — > oo) ~ x ( 2 - ff - A )/ A ^ an( j 

0(0, h) oc h (2 - a - &)/A = h l/s , with 6 = A/j3. (11) 

Finally, the isothermal susceptibility becomes 

§ = (S) =* ± |Tr,y = a + 2(A-l), (12) 

and upon eliminating A = /ft), one arrives at the follow- 
ing set of scaling relations 

a+f3{l +6) = 2 = or + 2/3 + r, T 1) • (13) 

Clearly, as consequence of the two-parameter scaling 
hypothesis Q, there can only be fwo independent ther- 
modynamic critical exponents. In the framework of 
Landau's mean-field approximation, the set of critical 
exponents reads a = 0, fj - \, y = 1,(5 = 3, and A = I; 
note that these integer or rational numbers really just 
follow from straightforward dimensional analysis. In 
both computer and real experiments, one typically mea- 
sures different critical exponent values, yet these still 
turn out to be universal in the sense that at least for 
short-range interaction forces they depend only on ba- 
sic symmetry properties of the order parameter and the 
spatial dimensionality d, but not on microscopic details 
such as lattice structure, nature and strength of inter- 
action potentials, etc. Indeed, the Ising ferromagnetic 
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and the liquid-gas critical points, both characterized by 
a scalar real order parameter, are governed by identical 
power laws, as is the critical behavior for planar mag- 
nets with a two-component vector order parameter and 
the normal- to superfiuid transition in helium 4, with a 
complex scalar or equivalently, a real two-component 
order parameter. The striking emergence of thermody- 
namic self-similarity in the vicinity of T c has been spec- 
tacularly demonstrated in the latter system, with the A- 
like shape of the specific heat curve appearing identical 
on milli- and micro-Kelvin temperature scales. 

2.3. Landau-Ginzburg-Wilson Hamiltonian 

In order to properly include the effects of fluctua- 
tions, we need to generalize the Landau expansion ([T|i to 
spatially varying order parameter configurations S(x), 
which leads us to the coarse-grained effective Landau- 
Ginzburg-Wilson (LGW) Hamiltonian 



WIS] 



d d x 



_s( X y + -[vs(x)Y 



+—S(xT-h(x)S(x) 



, (14) 



where r = a(T - T°) and u > as before, and h(x) 
now represents a local external field. Under the natu- 
ral assumption that spatial inhomogeneities are energet- 
ically unfavorable, the gradient term ~ [V5(x)] 2 comes 
with a positive coefficient that has been absorbed into 
the scalar order parameter field. Within the canonical 
framework of statistical mechanics, the probability den- 
sity for a configuration S (x) is given by the Boltzmann 
factor PAS] = exp(-<H[S]/k B T)/Z[h]. Here, the par- 
tition function Z[h] and expectation values of observ- 
ables A[S] are represented through functional integrals: 



ZW = §D[S]e-™l k » T , 
(A[S]} = J D[S]A[S(x)]T s [S] 



(15) 
(16) 



At h — 0, for example, the kth order parameter moments 
follow via functional derivatives 



k k 

'Y\S(-Xj)j = (k B T) k Y] 

' 7=1 



1 1 M -, Z H „• (17) 

1 f on{Xj) l/i=0 



and similarly the associated cumulants can be obtained 
from functional derivatives of ln^W; the partition 
function thus also serves as a generating function. For 
explicit calculations, one requires the integral measure 
in ( [To*) , e.g., through discretizing x — > x, on, say, a d- 
dimensional cubic hyperlattice, whence simply D[S ] - 



YlidS(xi). Alternatively, one may employ the Fourier 
transform S(x) = J S(q)e' q x ; noting that S(-q) 
S (q)* since S (x) is real, and consequently the real and 
imaginary parts of S(q) are not independent, one only 
needs to integrate over wave vector half-space, 



D[S] 



n 

<Ml>0 



d Re S (q) dlmS (q) 
V 



(18) 



In the Ginzburg-Landau approximation, one consid- 
ers only the most likely configuration S(x), which is 
readily found by the method of steepest descent for the 



path integrals in (16 1, leading to the classical field or 



Ginzburg-Landau equation 



= 



SW[S] 

6S(x) 



r- V 2 + -S(x) 2 
6 



S(x)-h(x). (19) 



In the spatially homogeneous case, ( 19 1 reduces to the 
mean-field equation of state Let us next expand in 
the fluctuations SS (x) — S (x) - <p about the mean order 
parameter <p — (S) and linearize, which yields 5h(x) m 
- V 2 + | <p 2 ^j 6S(x). Through Fourier transform one 
then immediately obtains the order parameter response 
function in the mean-field approximation, also known 
as Ornstein-Zernicke susceptibility 



xo(q) = 



dS(q) 



dh(q) 



1 



h=0 



1-2 



(20) 



where we have introduced the characteristic correlation 
length % = {r+l 4> 2 T l/1 , i.e., 

1/r l/2 

l/|2r| 1/2 



r>0 
r<0 



(21) 



On the other hand, consider the connected zero-field 
two-point correlation function (cumulant) 



C(x — x) 



(S(x)S(x'))-(S) 2 
5 2 \nZ[h] 



= (k B TY 



5h(x) 5h(x') 



(22) 



/!=() 



in a spatially translation-invariant system, we may de- 

j2^i C(q) e' q x , 

and through comparison with the definition of the sus- 
ceptibility in ( |20|) arrive at the fluctuation-response the- 
orem C(q) = k B T x{q), valid in thermal equilibrium. 
Generalizing the Ginzburg-Landau mean-field result 



(20 1, we may formulate the scaling hypothesis for the 



two-point correlation function in terms of the following 
scaling ansatz, which defines both the Fisher exponent r\ 
and the critical exponent v that describes the divergence 
of the correlation length £ at T c : 



C(T,q) = \q\- 2+ *C ± (qO, £ = £t M" 



(23) 
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The thermodynamic susceptibility then becomes 

X(t, q = 0) oc f-l cc\ T \- y , with y = v(2 -77) , (24) 

providing us with yet another scaling relation that con- 
nects the thermodynamic critical exponent y with rj and 
v. Consequently, we see that the thermodynamic critical 
point singularities are induced by the diverging spatial 
correlations. Fourier back-transform gives 



C(T,x) = \x\- (d - 2+ i ) C ± (x/OK{- 



(d-2+n) 



(25) 



at large distances \x\ — > co. In this limit, one expects 
(S (x) S (0)) — > <p 2 oc (-t) 2,8 , and comparison with ( 25 i 
therefore implies the hyperscaling relations 



(3 — - (d — 2 + rj) and 2 - a — dv . 



(26) 



The Ornstein-Zernicke function (|20b satisfies the seal 



ing law (|23|l with the mean-field values v = = and 



rj - 0. Notice that the set of mean-field critical expo- 



nents obeys (26 1 only in d — 4 dimensions. 



2.4. Gaussian approximation 

We now proceed to analyze the LGW Hamiltonian 



(14 1 in the Gaussian approximation, where non-linear 
fluctuation contributions are neglected. In the high- 
temperature phase, we have <p = 0, and thus simply omit 
the terms ~ uS(x) 4 , leaving the Gaussian Hamiltonian 



-I 



r + q'- 



\S(q)\ 2 -h(q)S(-q) 



(27) 



with the abbreviation f = J-^s- The associated Gaus- 
sian partition function is readily computed by complet- 



ing the square in ( 27 1, or the linear field transformation 
S(q) = S(q)-h(q)/(r + q 2 ), 



ZoW 



D[S] e 



-«o[S]/i B r 



I 

exp — — — — \Zo[h = 0], (28) 
\2k B T J q r + q l ) 



which yields the Gaussian two-point correlator 
(k B T) 2 (2n) 2d 6 2 Z [h] 



(S(q)S(q')) 



C (q)(2n) d 6(q + q'), C (q) 



ZoW 5h{-q)5h{-q') 
k B T 
r + q 2 



//=o 
■(29) 



Gaussian integrations give the free energy Fq[K\ = 
-k B T In ZoW of the model $2J}, 



Folk] 



2j q \r + 



+ k B TV In T 

r + q 1 



(30) 



Let us explore the leading singularity near in the spe- 
cific heat Ch=o - -T(d 2 Fo/dT 2 )i 1= o that originates from 
derivatives with respect to the control parameter r, 



C 



h=0 



V 



l<«vjr 



i 



(r + q 2 ) 2 ' 



(31) 



In high dimensions d > 4, the r.h.s. integral is 
UV-divergent, but can be regularized by a Brillouin 
zone boundary cutoff A ~ 27r/ao stemming from 
the original underlying lattice. Consequently the 
fluctuation contribution ( pT} is finite as r — > and 
a — as in mean-field theory. 

In low dimensions d < 4, we set k = q/ Vr = qij to 
render the fluctuation integral, which is UV-finite, 
dimensionless. With the c/-dimensional unit sphere 
surface area - 2jt d ^ 2 /F(d/2), one finds: 



C ft=0 k B (aT°) 2 ? 



■0\2 M-d r" 30 



f 

Jo 



iA-\ 



dk. (32) 



V ~ 2 d 7i d ' 2 T(d/2) Jo (l+k 2 ) 2 

As the critical temperature is approached, the cor- 
relation length prefactor diverges oc \T — T®\~~; 
already the lowest-order fluctuation contribution 
contains a strong infrared (IR) singularity that 
dominates over the mean-field power law. 

• At d — d c — 4, the integral diverges logarithmically 
as either A or ^ — > oo: 



Jo 



(l+k 2 ) 2 



dk ~ ln(A^) ; 



(33) 



note that at this upper critical dimension, ultravio- 
let and infrared divergences are intimately coupled. 

Above the critical dimension, we thus expect the mean- 
field scaling exponents to correctly describe the critical 
power laws. In dimensions d < d c , however, fluctuation 
contributions become prevalent and modify the mean- 
field scaling laws. Given the strongly fluctuating and 
correlated nature of critical systems, standard theoret- 
ical approaches such as perturbation, cluster, or high- 
and low-temperature expansions typically fail to yield 
reliable approximations. Fortunately, the renormaliza- 
tion group provides a powerful method to tackle inter- 
acting many-particle systems dominated by fluctuations 
and correlations, especially in a scale-invariant regime. 

2.5. Wilson 's momentum shell renormalization group 

The renormalization group program in statistical 
physics can be summarized as follows: The goal is to 
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establish a mathematical framework that can properly 
capture the infrared (IR) singularities appearing in ther- 
modynamic properties as well as correlation functions 
near a continuous phase transitions and in related situa- 
tions that are not perturbatively accessible. To this end, 
one exploits a fundamental new symmetry that emerges 
at a critical point, namely scale invariance, induced by 
the divergence of the dominant characteristic correla- 
tion length scale Approximation schemes need to 
carefully avoid the region where the physical IR sin- 
gularities become manifest; instead, one analyzes the 
theory in the ultraviolet (UV) regime, by means of ei- 
ther of various equivalent methods: In Wilson's mo- 
mentum shell RG approach, one integrates out short- 
wavelength modes; in the field-theoretic version of the 
RG, one explicitly renormalizes the UV divergences. 
Either method quantifies the weight of fluctuation con- 
tributions to certain coarse-grained or "renormalized" 
physical parameters and couplings. One then maps the 
resulting system back to the original theory given in 
terms of some "effective" Hamiltonian, which in Wil- 
son's scheme entails a rescaling of both control param- 
eters and field degrees of freedom. Thus one obtains re- 
cursion relations for effective, now scale-dependent run- 
ning couplings. Subject to a recursive sequence of such 
renormalization group transformations, these effective 
couplings will 

• either grow, and ultimately tend to infinity: to ac- 
cess a scale-invariant regime, one therefore has to 
set these relevant parameters to zero at the outset, 
which defines the critical surface of the problem; 

• or diminish, and eventually approach zero: these 
irrelevant couplings consequently do not affect the 
asymptotic critical scaling properties; 

• certain marginal parameters may also approach an 
infrared- stable fixed point, provided their initial 
value is located in the fixed point's basin of attrac- 
tion: clearly, scale-invariant behavior thus emerges 
near an IR-stable fixed point, and the independence 
from a wide range of initial conditions along with 
the automatic disappearance of the irrelevant cou- 
plings constitute the origin of universality. 

The central idea now is to take advantage of the 
emerging scale invariance at a critical fixed point as 
a means to infer the proper infrared scaling behavior 
from an at least approximative analysis of the ultravi- 
olet regime, where, e.g., perturbation theory is feasible. 
Thus one may establish a solid theoretical foundation 
for scaling laws such as |7]i and ( 23 1, thereby derive 



scaling relations, and also construct a systematic ap- 
proximation scheme to compute critical exponents and 
even scaling functions. We shall soon see that an appro- 
priate small parameter for a perturbational expansion is 
given through a dimensional expansion in terms of the 
deviation from the upper critical dimension e = d c — d. 

Wilson's momentum shell renormalization group ap- 
proach consists of two RG transformation steps: 

(1) Carry out the partition integral over all Fourier 
components S(q) with wave vectors residing in 
the spherical momentum shell A/b < \q\ < A, 
where b > 1: this effectively eliminates the short- 
wavelength modes. 

(2) Perform a scale transformation with the same scale 
parameter b > 1: x — » x' = x/b, q — > q' = bq. 
Accordingly, one also needs to rescale the fields: 



S(x) -> S'(x') = b ( S(x) , 
S(q)^S'(q') = bt- d S(q). 



(34) 



with a proper choice of ( ensuring that the rescaled 
residual Hamiltonian assumes the original form. 

Subsequent iterations of this procedure yield scale- 
dependent effective couplings, and the task will be to an- 
alyze their dependence on the scale parameter b. Notice 
the semi-group character of the above RG transforma- 
tions: there obviously exists no unique inverse, since the 
elimination step (1) discards detailed information about 
fluctuations in the UV regime. 

The mechanism and efficacy of the momentum shell 
RG is best illuminated by first considering the exactly 
tractable Gaussian model. Introducing the short-hand 
notations f < = f . .„ ^rL and f = f „ . . . 7 ^ s , 

Jq J\q\<A/b (2nf Jq JA/fc<| q|<A (2nf ' 

one may readily decompose the Hamiltonian into 
distinct additive Fourier mode contributions 



r + q^ 



\S(q)\ 2 -h(q)S(-q) 



Integrating out the momentum shell fluctuations then 
just gives a constant contribution to the free energy. 
We now wish to achieve that < Ho[S < ] — > 'H [S'] un- 
der the scale transformations in step (2). For the term 
q 2 \S(q)\ 2 , this is accomplished through the choice 



£ = ^ in (34 



ately result in t 



the other contributions then immedi- 
le following recursion relations for the 
r' = b 2 r, and h(q) — > 
-> h'(x') = b d ^h(x). 



control parameters r and h: r - 
h'(q') = b^h(q), whence h(x) 
Both the temperature variable r and the external field 
thus constitute relevant parameters, and the critical sur- 
face in parameter space is given by r — = h. As 



Uwe C. Timber / Nuclear Physics B Proceedings Supplement 00 (2011) 1 ^28\ 



7 



any other length scale, the correlation length scales ac- 
cording to £ — > f = £/b; eliminating the scale param- 
eter b one arrives at the relation far 



•1/2 



or v = s. 



Likewise, for the rescaled correlation function one finds 
C'(x') = £> 2 ^ C(x), whence 77 = 0: for the Gaussian the- 
ory, we recover the mean-field scaling exponents. 

We can gain additional non-trivial information by 
considering further couplings; e.g., imagine adding con- 
tributions of the form c s J d d x (y s S ) 2 to the Hamilto- 
nian ( 27 1 that represent higher-order terms in a gradient 
expansion for spatial order parameter fluctuations, sub- 
ject to preserving the Z2 and spatial inversion symme- 
tries. One readily confirms that c s — > c' s = b d ~ 2s ~ 2 ^c s = 
b~ 2( - s 'c s , which implies that all these additional cou- 
plings c s are irrelevant for s > 1 and scale to zero under 
repeated scale transformations. The inversion symme- 
try S(x) — > -S{x) permits general local non-linearities 
of the from u p J d d xS(x) 2p . Under Gaussian model RG 
transformation, these scale as u p — > u' p — b d ~ 2p ^u p = 
b 2p ~ (p ~ X)d u p ; these couplings are consequently marginal 
at d c (p) — 2p/(p — 1), relevant for d < d c (p), and ir- 
relevant for d > d c (p). The upper critical dimension 
decreases monotonously for p > 2, with the asymptote 
d c {oo) = 2. For the quartic coupling in the LGW Hamil- 
tonian ( fl4| i, this confirms d c (2) - 4, while d c (3) = 3 for 
a sixth-order term v J d d x S (x) 6 : v — > V = b 6 ~ 2d v, which 
becomes irrelevant near the upper critical dimension of 
the quartic term: v' = b~ 2 v at d c (2) — 4. In general, the 

coupling ratio -*P- = b 2 ~ d — renormalizes to zero in di- 

Up Up 

mensions d > 2. At two dimensions, the fields S(x) be- 
come dimensionless, £ = 0, and consequently all these 
non-linearities scale identically. The LGW Hamiltonian 
thus does not represent the correct asymptotic field the- 
ory, and one must resort to other effective descriptions 
(e.g., the non-linear sigma model). 

The above considerations already allow a discussion 
of the general structure of the momentum shell RG pro- 
cedure. According to (23 1, the general field rescaling 



(34|) should contain Fisher's exponent r\, £ = d ?" H? , 



whence h'(x) = U- d+2 ~ n ^' 2 h{oC), and (23 1 implies r' = 
b l l y T with the correlation length critical exponent v. In 
the simplest scenario, there are thus only two relevant 
parameters t and h. Let us further assume the presence 
of (a few) marginal perturbations — > u\ — u* + b~ x 'ui, 
while other couplings are irrelevant: v, — > v". = b~ y> Vi, 
with both Xj > and y t > 0. After a single RG transfor- 
mation step, the free energy density becomes 



/ S ing(T,/l,{M,},{v,})= (35) 

b- d ^ T y-^ + ^},{^}). 



After sufficiently many (£ » 1) RG transformations, the 
marginal couplings have reached their fixed point values 
m*, whereas the irrelevant perturbations have scaled to 
zero, 

fsm s (T,h,{Ui},{Vi}) = (36) 
Upon choosing b l |t| v = 1 for the scale parameter b e , one 



arrives at the scaling form 
fsms(r,K) = \r\ dv f ± (h/\i 



(v(d+2-Tj)l2 



(37) 



with f ± (x) = / s in g (±l,x, {m*},{0}). With the exponent 
identities ( [T3| and (26 1, this is equivalent to the scaling 
hypothesis (^|. In a similar manner, one readily derives 
the correlation function scaling law ( 25 1, employing the 
matching condition b e — £/g± for 

C(r, x, {«,}, {v,}) = b~ 2tl c(b clv T, {«*}, {0}) . (38) 

2.6. Dimensional expansion and critical exponents 

We are now ready to treat the non-linear fluctuation 
corrections by means of a systematic perturbation ex- 
pansion. The quartic contribution to the Hamiltonian 
reads in Fourier space 



74 



u f 

= * J*' 



S(qt)S (q 2 )S (qi)S (-qi -qi-qy) . (39) 

|<A 



Both the full partition function for the Hamiltonian ( 14 1 
and any associated TV-point correlation functions can 
then be rewritten in terms of expectation values in the 
Gaussian ensemble (we shall henceforth set kgT — 1) 

Z[h]=Zo[h](e-^ s A , 



v {UlSiq^e-^X 
\\S(qd)= ± r — r ^ . (40) 

Note that all expectation values of an odd number of 
fields S(qi) obviously vanish in the symmetric high- 
temperature phase, when the external field h = 0. Defin- 
ing the contraction of two fields as the Gaussian two- 
point function or propagator in Fourier space, 

S (q)S (q') = (S(q)S(q')} = C («) (2n) d 6(q + q') , 

we may write down Wick's theorem for Gaussian cor- 
relators containing an even number of fields, here a 
straightforward property of Gaussian integrations: 

(S(q l )S(q2)...S(q N - l )S(q N )}o = (41) 

/, S (<?ii(l))S (<?/ 2 (2)) ■ ■ ■ S (qi^^N-i^S (qi N (N)) ■ 



permutations 
i x {\)...i N (N) 
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Consequently, any arbitrary expectation value ( 40 1 can 
now be perturbatively evaluated via a series expansion 
with respect to the non-linear coupling u, and by means 
of (41 1 expressed through sums and integrals of prod- 



ucts of Gaussian propagators Co(g,). 

As an example, we consider the first-order fluctu- 
ation correction to the zero-field two-point function 
(S(q)S(q')) = C(q){2n) d 5{q + q'): (S(q)S(q')[l - % 
f lqil<A S(qi)S(q 2 )S(q 3 )S(-qi - q 2 - ?3)]} - According 
there are two types of contribu- 



to Wick's theorem (41 1 
tions: (i) Contractions of external legs S(q)S(q') yield 



terms that precisely cancel with the denominator in ( 40 1, 
leaving just the Gaussian propagator (S(q)S(q'))o. (ii) 
The twelve remaining contributions are all of the form 
f lqil<A S(q)S(qi) S(q2)S (q 3 ) S(-gi - q 2 - qi)S^ (q') = 



(2n) d 6(q + q')C (qfJ 
one obtains 



1pI<a 



Co(p). Collecting all terms, 



C(q) = C (q) 



u r 

■~C (q) C 
L ~Vl<A 



C (p)+O(u 2 ) 



(42) 



interpreting the bracket as the lowest-order contribution 
in Dyson's equation (see Chap.3.1 below), the integral 
turns out to be the associated self-energy to O(u), and 
d42l) can be recast in the form 



C{qT 



r + q + 



u f 

2 J\p\- 



1 



0(u 2 ) 



(43) 



)\p\<a r + p 1 - 

Notice that to order u, fluctuations here merely renor- 
malize the "mass" r, but there is no modification of 
the momentum dependence in the two-point correla- 
tion function C(q), implying that 77 will remain zero in 
this approximation. In a similar manner, one readily 
finds the first-order fluctuation correction to the four- 
point function at vanishing external wave vectors, i.e., 
the non-linear coupling u to be -\u 2 J^ <A Co(p) 2 . 

It is now a straightforward task to translate these per- 
turbation theory results into first-order recursion rela- 
tions for the couplings r and u by means of Wilson's 
RG procedure. To this end, we split the field variables 
in outer (5>: S{q) with A/b < \q\ < A) and inner (5<: 
S{q) with \q\ < A/b) momentum shell contributions; we 
then realize that there are four types of contributions: 

• terms involving merely inner shell fields that are 
not integrated, e.g. ~ u J S A < e~' Ho ^ just need to 
be re-exponentiated; 

• integrals such as u ^ S\S > e^ ols ^ vanish; 

• contributions ~ u J > S>e~' Ho ^ that contain only 
outer shell fields become constants that directly 
contribute to the free energy; 



for terms of the form ~ u f S< 5>e _Wo , one has to 
perform Gaussian integrations over the outer shell 
fields 5>, yielding corrections to the propagator for 
the inner shell modes. 



Employing ( [43J >, using 77 = 0, and introducing Sd 
K d /(2n) d = l/2 d - 1 ^ /2 r(c//2), one thus finds to 0(u): 



r'=b 2 



Mr) 



r + 



b 4 - d u 



2 J A/b r + 



dp 



3m 



A/b 

B(r) 



= b 



A-d, 



L J A 



A pd -X 



dp 



A/b (r + p 2 ) 2 



(44) 



(45) 



For T » T c , or r — > 00, the fluctuation corrections be- 
come suppressed, and one recovers the recursion rela- 



tions r' = b r and u' = b u of the Gaussian theory. 
Near the critical point, i.e. for r« 1, one may expand 

1 - b 2 - d 



A(r) = S d A 



d-2 



d-2 



-rS d A d ' A 



B(r) = S d A 



d-4 



1 



d-4 

b 4-d 



d-4 



+ 0(r 2 ) , (46) 
+ O(r) . (47) 



It is useful to consider instead differential RG flow 
equations that result from infinitesimal RG transforma- 
tions. Setting b = e se with 6( -> 0, (M-(|47} turn into 



— _2r(() + —S d A 



f(()u({) 



du({) 
dl 



(4 - d) u{{) 
+0(ur, ii 2 ) . 



S d A d - 4 + 0(ur 2 ,u 2 ) , (48) 
3 



u({) 2 S d A d - 4 



(49) 



We specifically seek renormalization group fixed points 
(r*,u*) that describe scale-invariant behavior, to be de- 
termined by the conditions dr(()/d( — = du{€)jdl. 
There is obviously always the Gaussian fixed point 
= 0; linearizing ([49]) in terms of the deviation 



6uo({) = u({)-Uq, one finds dduo({)/d{ « (d — 4)5uq{€); 
u* Q is hence stable for d > d c — 4, but unstable for d < 4. 
Below the upper critical dimension, there exists also a 
positive Ising fixed point u*S d = | (4 - d)A A ~ d , which 
is then also stable since d6iij({)/d{ ~ (4 - d)6iii{€) for 
6u\(l) = u{€)—u\. Correspondingly, the critical behavior 
is governed by the Gaussian fixed point and associated 
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scaling exponents in dimensions d > 4, but by the non- 
trivial Ising fixed point in low dimensions d < d c = 4. 
Notice also that the numerical value of the Ising fixed 
point becomes small near the upper critical dimension, 
and indeed u* emerges continuously from ul — as 
e = 4 - d is increased from zero. At the non-trivial 
RG fixed point, e may serve to provide a small effective 
expansion parameter for the perturbation expansion. 

To lowest order in the non-linear coupling, ( |48] i yields 
at the Ising fixed point r* = - \ u\S dA d " 2 = -\ eh 2 , 
which describes a non-universal, fluctuation-induced 
downward shift of the critical temperature. Introducing 
the deviation t = r-r\ - a(T - T c ) from the true critical 
temperature T c , one may rewrite the flow equation (48 i 
to obtain the recursion relation for this modified relevant 
running coupling: 



df(€) 
dl 



2 



(50) 



Its solution in the vicinity of the Ising fixed point reads 
?(£) = f(0)exp^2 - |^J. Combining this result with 

£,(€) = £(0) e~ e , one identifies the correlation length ex- 
ponent v = 2 - |, or, in a consistent expansion to first 
order in e = 4 - d: 

v=\ + -^ + 0(e 2 ), 77 = + <9(e 2 ). (51) 

As anticipated, the critical exponents depend only on 
the spatial dimension, not on the strength of the non- 
linear coupling or other non-universal parameters. Note 



that at d c = 4, (|49j is solved by u{€) - ft(0)/[l + 
3u(Q)£/16n 2 ], a very slow approach to the Gaussian 
fixed point that induces logarithmic corrections to the 
mean-field critical exponents, see ( |8T| ) below. We finally 
remark that the RG procedure generates novel coupling 



terms ~ S , V S , etc. To order e , their feedback into 
the recursion relations can however be safely neglected. 

In summary, the renormalization group procedure as 
outlined above in Wilson's momentum shell formula- 
tion allows us to derive hitherto phenomenological scal- 
ing laws, and thereby gain deeper insights into scale- 
invariant features. We have also seen that the number 
of relevant couplings (two at standard critical points, 
namely t and h) equals the number of independent crit- 
ical exponents. Below the upper critical dimension, 
fluctuation corrections modify the critical scaling dras- 
tically as compared to the mean-field predictions. Fi- 
nally, perturbative calculations that are safely carried 
out in the UV regime may be employed to systemati- 
cally compute scaling exponents through a power series 
in the dimensional parameter e — d c — d. 



3. Field Theory Approach to Critical Phenomena 

While Wilson's momentum shell scheme renders the 
basic philosophy of the renormalization group transpar- 
ent, it becomes computationally quite cumbersome once 
nested momentum integrals appear beyond the first or- 
der in perturbation theory. Unnecessary technical com- 
plications in evaluating fluctuation loops can be avoided 
by extending the UV cutoff to A — > oo, at the price of 
divergences in dimensions d > d c . However, we already 
know that the Gaussian theory governs the infrared 
properties in that dimensional regime; thus in statisti- 
cal physics theses UV singularities do not really pose a 
troublesome issue. We may however employ powerful 
various tools from quantum field theory, proceed to for- 
mally renormalize the UV divergences, and thereby gain 
crucial information about the desired IR scaling limit, 
provided an IR-stable RG fixed point can be identified 
that allows us to connect the UV and IR regimes. This 
chapter provides a succinct overview of how to con- 
struct the perturbation expansion in terms of Feynman 
diagrams for one-particle irreducible vertex functions, 
proceeds to analyze the resulting UV singularities, and 
finally utilizes the renormalization group equation to 
identify fixed points and determine the accompanying 
critical exponents. For more extensive treatments of the 
field-theoretic RG approach to critical phenomena, see 
Refs. ||T0l-|[l5l and other excellent texts. 

3.1. Perturbation expansion and Feynman diagrams 

We now generalize our analysis to a LGW Hamilto- 
nian with continuous O(n) order parameter symmetry 



+ 



S a (x) 2 



1 



[VS a (x)] 2 



(52) 



which encapsulates the critical behavior for the Heisen- 
berg model for a three-component vector order param- 
eter (n = 3), the planar XY model (and equivalently, 
superfluids with complex scalar order parameter) for 
n = 2, reduces to Ising Z2 symmetry for n = 1, and 
in fact describes the scaling properties of self-avoiding 
polymers in the limit n — > 0. As in (40i, one constructs 



the perturbation expansion for arbitrary A-point func- 
tions (FJi m terms of averages within the Gaussian 
ensemble with u — (keeping k^T = 1). Diagrammat- 
ically, the Gaussian two-point functions or propagators 
Co(q) 5 a P = 6 a P l(r+q 2 ), diagonal in the field component 
indices, are represented through lines, to be connected 
through the non-linear vertices — | : 
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In the presence of external fields h", the partition 
function serves as generating functional for correlation 
functions (cumulants): 

Z[h] = ^exp Jd d xY J h a S a 

n°%=n s ™ zlh] 



5h a - 



\h=0 



(53) 



The cumulants are graphically represented through con- 
nected Feynman diagrams; e.g., for the propagator: 



Note that the second graph is a mere repetition of the 
first; indeed, upon defining the self-energy E as the sum 
of all one-particle irreducible Feynman graphs that can- 
not be split into lower-order contributions simply by 
cutting a propagator line, one infers the following gen- 
eral structure for the full propagator C(q): 

— + •■• 



The second line is a graphical depiction of Dyson's 
equation that reads in Fourier space C(q) - C[)(q) + 
C (q)Hq)C(q), solved by C{qT x = C {qT X -2(g). 

In order to similarly eliminate redundancies for arbi- 
trary Appoint functions, one proceeds with a Legendre 
transformation to construct the generating functional for 
vertex functions: 



n*] 



= - In Z[h] + fd d x h a a> a , 



n 



(JO"' \h=a 



(54) 



where O ff = 6lnZ[h]/6h a . Through appropriate func- 
tional derivatives, these vertex functions can be related 
to the corresponding cumulants, for example for the 
two- and four-point functions: 

T^ 2 \q) = C{qT X , 



Y\S(qd) =-flcfe)r (4) (te}). (55) 



By means of these relations one easily confirms that 
the perturbation series for the vertex functions precisely 



consist of the one-particle irreducible Feynman graphs 
for the associated cumulants. Moreover, the perturba- 
tive expansion with respect to the non-linear coupling 
u diagrammatically translates to a loop expansion. Ex- 
plicitly, one obtains for the two-point vertex function to 
two-loop order (resulting from the first, third, and fourth 
graph above): 



I< 2 \q) = r + q 2 + 



n + 2 



2 f 1 

«f f —1 f — 

/ Jkr + k 2 J k , ( r + k' 



(56) 



n + 2 
18 



M 2 r ' r _j i 

U Jkr + k 2 Jtr + k' 2 r + (q-k- k') 2 ' 



for the four-point vertex function to one-loop order, we 
just have the single Feynman diagram 



which yields at vanishing external wave vectors 

r< 4 >( te = o}) = «- n -^u 2 f (57) 

3.2. UV and IR divergences, renormalization 

Let us now investigate the fluctuation correction ( f57| 
to the four-point vertex function. In dimensions d < 4, 
we can safely set the UV cutoff A — > oo, and obtain after 
rendering the integral dimensionless: 



r d d i 

u I 

J (2n 



1 



) d (r + k 2 ) 2 

Z+d/2 



2««-V/ 2 r(d/2) Jo (l+x 2 ) 2 



f 

Jo 



■ dx . 



(58) 



Observe that the true effective coupling in the perturba- 



tion expansion is not u, but u r' 



(d-4)/2 



oo as r 



0, 



indicating the emerging infrared divergences that ren- 
der a direct perturbative approach meaningless in the 
critical regime. We conclude once again that for d < 
d c = 4, fluctuation corrections are IR-singular, and con- 
sequently expect the critical power laws to be modified 
as compared to the mean-field or Gaussian approxima- 
tions. In contrast, the integral remains regular in the 
infrared for d > 4, but becomes ultraviolet-divergent. 
Indeed, keeping the cutoff finite, we have in dimensions 
larger or equal to the upper critical dimension d c — 4: 



X 



ud-\ 



o (r + k 2 ) 2 



dk 



ln(A 2 /r) 



d = A 
d>4 



(59) 



which both diverge as A — > oo. Notice that at d c — 
4, the logarithmic IR and UV divergences are coupled, 
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signaling the scale invariance of the LGW Hamiltonian 
in four dimensions. 

One may take a convenient shortcut to determine 
the critical dimension through simple power counting 
in terms of an arbitrary momentum scale p. Lengths 
then scale as [x] = pT 1 , wave vectors as [q] = p, and 
fields consequently have the (naive) scaling dimension 
[S a (x)] = = p- l+dl2 . As noted before, the fields be- 
come dimensionless in two dimensions. For the LGW 



Hamiltonian ( 52 1 with continuous rotational symmetry 
in order parameter space, that is also just the lower 
critical dimension: For d < d\ c — 2, the system can- 
not maintain spatially homogeneous long-range order 
if n > 3; and at die = 2 and for n — 2 only quasi- 
long-range order may exist with algebraically decay- 
ing correlations and temperature-dependent decay ex- 
ponent (Berezinskii-Kosterlitz-Thouless scenario). For 
the couplings in (52 1 one infers [r] = p 2 , which means 
that the temperature control parameter constitutes a rel- 
evant coupling, and [u] = p A ~ a ', which is relevant for 
d < 4, marginal at d c = 4, and (dangerously) irrelevant 
for d > 4. Furthermore, dimensional analysis confirms 
that fluctuation loops become UV-divergent only for the 
vertex functions r (2) (g), but only up to order q 2 in a 
long-wavelength expansion, and F (4) ({g, = 0}). The fol- 
lowing table summarizes the mathematical distinctions 
and their physical implications in the different dimen- 
sional regimes. 



dim. 


perturb. 


0(n) (D 4 


critical 


range 


series 


field theory 


behavior 


d<2 


IR-sing. 


ill-defined 


no long-range 




UV-conv. 


u relevant 


order (n > 2) 


d<4 


IR-sing. 


super-ren. 


non-classical 




UV-conv. 


u relevant 


exponents 


d c =4 


log. IR-/ 


renorm. 


logarithmic 




UV-div. 


u marginal 


corrections 


d>4 


IR-reg. 


non-renorm. 


mean-field 




UV-div. 


u irrelevant 


exponents 



It is useful to perform the loop integrations in dimen- 
sional regularization; i.e., to assign the following values 
to wave vector integrals, even for non-integer d and <x: 



! 



d d k k 2ir 

(liff (t + k 2 ) s ~ 
r(o- + d/2)T(s-cr-d/2) 
2 d n d / 2 T(d/2)T(s) 



(60) 



JT-S+djl 



For d < d c , where no UV divergences appear, this result 
follows directly by introducing spherical coordinates in 
momentum space. Beyond the upper critical dimension, 



essentially divergent surface integrals are discarded in 
( |60] >. UV singularities become manifest as dimensional 
poles in Euler's F functions. 

We may now proceed with the renormalization pro- 
gram. The goal is to absorb the UV divergences into 
renormalized couplings that through this procedure be- 
come scale-dependent. Beginning with the order pa- 
rameter susceptibility, we naturally demand that;^ -1 = 
C(q = 0)~ l = T®\q = 0) = t = r - r c , i.e., x diverges 
at the true critical temperature T c . From (56 1 we thus 
obtain to first order in u, 



n + 2 
6 

n + 2 



u I 

Jk 



1 



uS, 



kr c +k 2 

A d-2 



+ 0(u 2 ) 



d-2 



+ 0(u 2 ) , 



(61) 



which is to be interpreted as a non-universal fluctuation- 
induced downward shift of the critical temperature, the 
analog of r* in Wilson's scheme. As A — > oo, r c 
becomes quadratically UV-divergent near four dimen- 
sions; this divergence becomes absorbed into the new 
proper temperature variable t by means of an additive 



renormalization. Inserting (61 1 into (56 1 then yields to 
order u 



x(qT 



q + T 



1 



n + 2 C 



i 



k 2 (T + k 2 ) 



(62) 



At the upper critical dimension, ( |62| ) and §57) are log- 
arithmically divergent as A — > oo, showing up as 1/e 
poles in the dimensionally regularized integral values. 
These UV poles are subsequently absorbed into renor- 
malized fields S'x and parameters through the following 
multiplicative renormalization prescription: 



Z l J 2 S a 



T R = Z r Tp 



-N/2 



u R = Z u uA d p' 



d-A 



(63) 
(64) 



which defines dimensionless renormalized couplings t r 
and u R , and where = F{3-d/2)/2 d ~ l tv 4 ^ 2 is a regular 
(near d c = 4) geometric factor. To this end, one must 
carefully avoid the IR-singular regime, i.e., evaluate the 
fluctuation integrals at a safe normalization point, e.g., 
tr — I (or q = p). In the minimal subtraction procedure, 
the renormalization constants Z s in ( |63~| and Z T , Z u in 
( |64") i are chosen to contain only the 1/e poles and their 
residua. This leads to the following Z factors 



0(u R ) : Zr = 1 



Z, = 1 - 



n + 2 Ur 
6 T 

n + 8 Ur 



0{u\) : Z s = 1 + 



n + 2 u 



144 



(65) 
(66) 

(67) 



Uwe C. Tauber I Nuclear Physics B Proceedings Supplement 00 (2011) I- \28\ 



12 



all calculated to first non-trivial order in u R by means 



of dimensional regularization ( 60 1 and within the mini- 
mal subtraction prescription. Z T and Z„ follow directly 
from the one-loop results ( 62 1 and ( 57 1, whereas Z$ = 1 



to 0(u R ) due to the absence of any wave vector de- 
pendence in the "Hartree" loop, whence field renormal- 
ization only ensues to two-loop order from the right- 
most "sunset" Feynman diagram in the propagator self- 



energy, and the 1/e pole in the final expression in (56 1 



The first-order linear partial differential equation ( 69 1 
can next be formally solved via the standard method of 
characteristics; to this end, one lets p — > p{t) — p C, 
with a dimensionless scale parameter i ; note that in con- 
trast to the convention in Chap. 2.5, the IR regime is now 
reached in the limit I — > 0. Inserting this parametriza- 



tion into the RG equation (69 1, one obtains an equivalent 
set of coupled first-order ordinary differential equations, 
namely the RGflow equations for the running couplings 



3.3. RG equation and critical exponents 

Through the selection of a normalization point well 
outside the critical regime, the renormalized fields and 



parameters in ( 63 I, ( 64 1 explicitly carry the momentum 



scale p.. On the other hand, the unrenormalized quan- 
tities, including the Appoint vertex functions, naturally 
do not depend on this arbitrary scale p: 

= ± TW(t, u) = L [ Z f 2 Tf (p, TR , u R )] . (68) 



dp 



Carrying out the derivative with respect to p by taking 
into account the scale dependence of Z s , t r , and u R , (68 i 
can be rewritten as a partial differential equation, 



d 



a 



d 



N 

H + 7T 7s + Jt t r — +Pu -5— 
op 2 ot r ou R 

lf ) ip,T R ,U R ) = 0. 



(69) 



This Gell-Mann-Low renormalization group equation 
carries crucial information on the fundamental scale de- 
pendence of the renormalized physical system, here ren- 
dered explicit for the vertex functions. In ( |69| we have 
introduced Wilson's flow functions defined as 



Ys 



d \ 1 7 
: ^lo lnZ * 



n + 2 
72 



u\ + 0{u\) 



d I , t r 
■ p — In — 

^ dp lo T 



n + 2 



u R + 0(u R ) 



(70) 



(71) 



where the second lines follow from the lowest-order re- 
sults ( |6*7| i and ( |65) , and the RG beta function for the 
non-linear coupling u, 



d I 



u R 



d-4 + 



lnZ " 



= UR 



-e + 



u R + 0(u 2 R ) 



(72) 



where d66|l has been inserted. 



di dl 



= Pu(V, (73) 



with initial values f(l) = r s , S(l) = u R , and simi- 
lar differential equations for the Appoint vertex func- 
tions, which involve their overall naive scaling dimen- 
sions and the anomalous contributions stemming from 
the field renormalization, as encoded in d70j. 

For example, for the susceptibility x(l) = r' 2 '^)" 1 , 
onehas^ R (^,r fi ,M S ,^) _I = p 2 x R (j R ,u R ,q/py\ and its 
RG flow correspondingly integrates to 



XR((r l =x R ar l e 2 ex P 



r 



di' 



(74) 



Near an infrared-stable RG fixed point u*, i.e., a zero 
of the RG beta function fi u (u*) = with f3' u (u*) > 0, 
the flow equation for the running temperature variable 
is readily solved: ?({) « t r (?" t , where y* = y T (u*). In- 
serting this and the fixed point value y* s = ys(u*) into 
(|74"1> yields the following general scaling form 



'p 2 C 2+ ^MrR^,u,qlper x 



(75) 



With the matching condition E = \q\lp, one thus recov- 



ers (23 1, with the identifications rj = -y* and v = -l/y* 
for the critical exponents. 

Our perturbative RG analysis of the (9(n)-symmetric 



LGW Hamiltonian (52 1 yielded the one-loop beta func- 



tion ( 72 1, whose zeros are (i) the Gaussian fixed point 



0, stable for e < or d > d c = 4, obviously result- 



ing in the mean-field exponents 77 = and v 
(ii) the non-trivial Heisenberg fixed point 



6e 
« + 8 



+ 0(e z ) 



and 



(76) 



which exists and becomes IR-stable for e > 0, i.e., in 
dimensions d < 4. This allows us to compute the critical 
exponents in a systematic e = 4 — d expansion, 



n + 2 
2(n + 8) ; 



e 2 + 0(e 3 ) 



n + 8 



(77) 
(78) 
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Aside from the dimensionality, these values only de- 
pend on the number of order parameter components n. 
Note that ( |78] > reduces to the Ising exponent values ( |5T) 
for n = 1; in the limit n — > oo one also finds the cor- 
rect exponents for the exactly solvable spherical model, 
namely rj — and v = 1/(2 - e) = l/(d - 2), which 
diverges at the lower critical dimension d\ c = 2. At the 
upper critical dimension d c = 4, the solution of the flow 
equation for the running non-linear coupling reads 



u{€) = 



"R 



whence approximately 
tr 



(79) 



£ 2 (\n |^|)(«+2)/(»+8) ' 



(80) 



which in turn implies that the correlation length diver- 
gence picks up logarithmic corrections to the mean-field 
behavior, 



^r s I/2 (lnT*) 



(n+2)/2(n+8) 



(81) 



The field-theoretic formulation of the renormaliza- 
tion group provides us with an elegant and powerful 
tool to extract the proper infrared scaling properties in 
low dimensions d < d c from a continuum theory via 
a careful analysis of its ultraviolet singularities that ap- 
pear in dimensions d > d c . The renormalization group 
equation carries information on the scale dependence of 
physical parameters and correlation functions, and al- 
lows to make connections between the UV and IR lim- 
its provided a stable RG fixed point can be identified. 
One may then systematically derive scaling laws, and 
acquire a thorough understanding of the origin of uni- 
versality and its realm of validity for a specific phys- 
ical system. Moreover, we have seen that a perturba- 
tive analysis allows a controlled computation of critical 
exponents (and also of scaling functions) in the frame- 
work of a dimensional expansion near the upper criti- 
cal dimension. This e expansion certainly provides use- 
ful information on overall trends, and can in some in- 
stances even be rendered to a precision calculation if 
sufficiently high orders in the perturbation series can be 
evaluated and subsequently be refined through Borel re- 
summations. In addition, modern non-perturbative "ex- 
act" numerical RG methods have succeeded in yielding 
very accurate results. It should also be stressed that it is 
of course the very concept of universality that also al- 
lows us to infer meaningful information from numerical 
simulations of simplified lattice or continuum models. 



4. Critical Dynamics 

We now venture to investigate dynamic critical phe- 
nomena near continuous phase transitions, first in the 
vicinity of critical point in thermal equilibrium, and 
later at genuine non-equilibrium phase transitions in ex- 
ternally driven, non-isolated systems. The natural time 
scale separation between the slow kinetics of the or- 
der parameter (along with any other conserved fields) 
and fast, non-critical degrees of freedom suggests a 
phenomenological description in terms of non-linear 
stochastic Langevin-type differential equations, and al- 
lows a generalization of universal scaling laws to time- 
dependent phenomena. Distinct dynamical universal- 
ity classes ensue dependent on the order parameter it- 
self representing a conserved quantity or not, and po- 
tentially its dynamical coupling to other conserved hy- 
drodynamic modes [16|. We begin by writing down 
scaling laws for dynamical response and correlation 
functions, and subsequently introduce effective meso- 
scopic Langevin equations, with stochastic noise cor- 
relations that near thermal equilibrium are constrained 
by fluctuation-dissipation relations. It is then demon- 
strated how such non-linear stochastic partial differen- 
tial equations can be mapped onto a field theory repre- 
sentation via the Janssen-De Dominicis response func- 
tional, which in turn may be analyzed by means of the 
field-theoretic renormalization group tools developed in 
the previous chapter. We will specifically construct a 
dynamic perturbation theory expansion and determine 
the universal scaling behavior along with the dynamic 
critical exponents for the relaxational models A and B 
with non-conserved and conserved order parameter, re- 
spectively; for more in-depth treatments, the reader is 
referred to Refs. ED-lDI) and chaps. 4,5]. In addi- 
tion to purely relaxational kinetics, we shall also address 
the critical dynamics of isotropic ferromagnets GUI , as 
well as generic scale invariance and non-equilibrium 
phase transitions in driven diffusive systems ETI . 

4.1. Dynamical scaling hypothesis 

Let us first recall the behavior of the static order pa- 
rameter correlation function and susceptibility near a 
critical point (h = and r — > 0), captured by the scaling 



laws ( 23 i and ( 25 i, and induced by the divergence of the 
characteristic correlation length, ^(r) ~ \t\~ v . As spa- 
tially correlated regions grow tremendously upon ap- 
proaching the phase transition, one expects the typical 
relaxation time associated with the order parameter ki- 
netics to increase as well, t c (j) ~ %(j) z ~ |r|~ zv . This 
phenomenon of critical slowing-down is governed by 
the dynamic critical exponent z — v,/v, which can also 
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be visualized as the ratio between the exponents that de- 
scribe the divergence of correlations in the temporal and 
spatial "directions", respectively. We may thus write 
down a dynamic scaling ansatz for the corresponding 
wavevector-dependent characteristic frequency scale, 



io c (j,q) = \qf&±{q%) 



(82) 



with co±(y — > oo) — > const., whence the critical disper- 
sion relation becomes co c (0,q) ~ \q\ : . 

We are particularly interested in describing the time 
dependence for the order parameter response and corre- 
lation functions: 



X(x-x',t-t') = 



d(S(x,t)) 



dh(x',t') 



h=0 



C(x, t) = (S(x,t)S (0,0)> -(Sf 



(83) 
(84) 



where we are considering a stationary dynamical regime 
where spatial and temporal time translation invariance 
holds. In thermal equilibrium (only !), the spatio- 
temporal Fourier transforms of these functions are in- 
timately connected through the fluctuation-dissipation 
theorem (FDT) 



C(q, a>) = 2ks T Im 



x(q> w) 



to 



(85) 



Generalizing the static scaling laws ( |23) and ( |25] >, we 
may then formulate the dynamical scaling hypothe- 
sis for the asymptotic critical properties of the time- 
dependent susceptibility and correlation function: 



C(t, x, t) = \x\- (d - 2+,,) C ± (x/£, t/t c ) . 



(86) 
(87) 



As a consequence of the stringent constraints imposed 
by the FDT ( |85j ), the same three independent critical 
exponents v, 77, and z characterize the universal scaling 



regimes in both (86 1 and (87 1. Away from thermal equi- 
librium, where the FDT restrictions do not apply, the dy- 
namic response and correlation functions are however 
in general characterized by distinct scaling exponents. 
We remark that appropriate variants of the dynamical 
scaling hypothesis may also be postulated for transport 
coefficients. 

4.2. Langevin dynamics and Gaussian theory 

The critical slowing-down of the order parameter ki- 
netics produces an effective time-scale separation be- 
tween the critical degrees of freedom, additional con- 
served hydrodynamic modes that might be present, and 
all other comparatively "fast" variables. This obser- 
vation naturally calls for a mesoscopic Langevin de- 
scription of critical dynamics, where the fast degrees 



of freedom are treated as mere white noise that ran- 
domly affects the few slow modes in the system. In 
such a coarse-grained picture, one writes down cou- 
pled stochastic equations of motion for the order pa- 
rameter and perhaps any other conserved fields that re- 
flect their intrinsic microscopic reversible dynamics as 
well as irreversible relaxation kinetics, the latter con- 
nected in thermal equilibrium with the noise strengths 
through Einstein relations or FDTs. Generally the var- 
ious possible mode couplings of the order parameter to 
additional conserved, and consequently diffusively slow 
modes leads to a splitting of the static into several dy- 
namic universality classes ifTBI [TBI |9"1 . 

Here we shall assume that the order parameter field 
is decoupled from any other slow modes, and first fo- 
cus on its purely relaxational kinetics |9] Chaps. 4,5]. 
If the order parameter itself is not a conserved quantity, 
any deviation from thermal equilibrium will just tend to 
relax back to the minimizing configuration of the free 
energy, e.g. given by the (9(n)-symmetric LGW Hamil- 
tonian 1 



(88) 



Bt 



5S a (x,t) 



with Gaussian white noise that is fully characterized by 
its first two moments, 



(£*{x,t)if(x',t')) = 

2Dk B T6(x - x')5{t - t')6 al3 



(89) 



As can be inferred from the associated Fokker-Planck 
equation, the Einstein relation that connects the noise 
correlator strength with the relaxation constant D and 
temperature guarantees that the probability distribu- 
tion for the field S a asymptotically approaches the 
canonical stationary distribution , t] — > !Pj[5*] °c 
exp(-'7Y[5 > ]/^ B r) as t -> 00. 

If the order parameter is conserved under the dynam- 
ics, satisfying a continuity equation, its spatial fluctua- 
tions can only relax diffusively; as a consequence, one 
needs to replace the constant relaxation rate D by the 
diffusion operator -DV 2 , both in the Langevin equa- 
tion ([88) and the noise correlation ( |90) . In the follow- 
ing, we shall treat both these situations simultaneously, 
letting D — > D(iV) a , where either a = 0, corresponding 
to the purely relaxational model A of critical dynamics 
as appropriate for a non-conserved order parameter; or 
a = 2, which describes model B with a conserved order 
parameter field. 

Let us again begin with the Gaussian or mean-field 
approximation, where we set the non-linear coupling 
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u — 0. A Fourier transform in space and time accord- 
ing to S a (x, t) = f-^sj ^S(q,oS)e^ x - iM of the thus 
linearized Langevin equation (|88]l yields 



[-KJ + D^ fl (r + ? 2 )] S a (q, co) = 

Dq a h a (q, co) + f(q, to) , 

with {C"{q, co)) = and 

(t a (q,co)f(q',co')) = 2kBTDq a 

(2n) d+y 5(q + q')5(oj + oj')8 al3 



(90) 



(91) 



and where the external field term -J] a h a S a has been 
added to the LGW Hamiltonian ( |52) , A derivative of 
( 90 1 with respect to the external field then immediately 
gives the dynamic response function 



ap d(S a (q,co)) 
X ( f(q, to) - 



G (q, co) = 



dhP{q,co) Ia=o 

Dq a G (q,to)6 a P , 
1 

-ico + Dq a (r + q 2 ) 



(92) 



Its temporal Fourier backtransform of course satisfies 
causality: Go(q, f) vanishes for t < 0, and reads 



G (q, t) = e 



-Dq\r+q 2 )t 



0(f) 



(93) 



from which we infer the characteristic relaxation rate 
t~ l = Dq a (r + q 2 ): For model A (a - 0), the order 
parameter relaxes diffusively at the critical point (i.e., 
z = 2), while for model B {a = 2) critical slowing-down 
induces a crossover from Drq 2 to Dq A as r — > (z = 4). 
For h a = 0, the dynamic correlation function is readily 
obtained from |90| and ( |9"T| i: 

(S a (q,co)S fi (q',co')) = C (q,co) 
(2n) d+l 5{q + q')S(co + co')6 a0 , 
2k B TDq a 



C (q, co) 



C (q, t) 



co 2 + [Dq a (r + q 2 )] 2 
= 2k B TDq a \G Q (q,C0)\ 2 , 

_ k*r_ 

r + q 2 



-Dq° (r+q- 



(94) 



(95) 



Comparing these results with ([86j» and ( |87j i, one again 
identifies the static Gaussian critical exponents v — \ 
and r\ — 0, and the mean-field dynamic exponent z — 
2 + a for the purely relaxational models A and B. 

4.3. Field theory representations of Langevin dynamics 

This subsection describes how stochastic Langevin 
equations of motion can be mapped onto continuous 
field theory representations. To this end, we consider 



the following general coupled Langevin equations for 
some mesoscopic stochastic variables S a : 



dS a (x,t) 



= F a [S](x,t) + r(x,t), (96) 

(97) 



dt 

(C(x,t)cf(x',t')) = 

2L a [S]S(x-x')S(t-t')6 afl 



Naturally we assume (£ a (x, t)) = here, since a non- 
vanishing mean of the stochastic forces or noise could 
just be included in the systematic forces F a [S]. Note 
that the noise correlator L a may be an operator, as is 
the case for conserved variables, and could also be a 
functional of the slow fields S a . The crucial input is 
that we assume the noise history to represent a Gaus- 
sian stochastic process, whose probability distribution 
if completely determined by the second moment (|97|): 



oc exp 



1 -fd d x£dtY J r(L a y 1 F 



(98) 



Switching dynamical variables from the noise g a to 
the slow stochastic fields S a yields < W[£\D[£] = 
P[S]D[S] oc e~ S[S] D[S], with the Onsager-Machlup 
functional providing the associated exponential weight 
that may be viewed as a field theory action: 

g[S] = l -J d d x£dtJ^(d,s a -F a [S]) 

[(L a y l (d,S a -F a [S])]. (99) 

The observant reader will have noticed that the func- 
tional determinant stemming from the variable transfor- 
mation has been ignored here; however, upon utilizing 
a forward time discretization, i.e., the Ito interpretation 
for non-linear stochastic processes, this functional de- 
terminant turns out to be constant, and can simply be 
absorbed into the functional integral measure. Notice 
also that the overall nor malization J £>[£]W[£] = 1 im- 
plies the corresponding "partition function" to be unity, 
and hence to carry no information, in stark contrast 
with thermal equilibrium statistical mechanics. While 



the Onsager-Machlup functional ( 99 1 provides a desir- 
able field theory representation of stochastic Langevin 
dynamics, it is also plagued by two technical prob- 
lems: First, it contains (L a )~ l , which for conserved 
fields entails an inverse differential operator or Lapla- 
cian Green's function; second, it includes the square 
of the systematic force terms F Q '[5] and consequently 
highly non-linear contributions. It is thus beneficial 
to partially linearize the above action by means of a 
Hubbard-Stratonovich transformation, at the expense of 
introducing an additional dynamical field variable. 
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In order to completely avoid any possible singular- 
ities incorporated in the inverse operator (L*) , we 
follow another more direct route here. The goal is 
to compute averages of observables A that should be 
functionals of the slow modes S a over noise "histo- 
ries": (A[S]} ( oc fD[£]A[S(fj]W[£]- Inserting a 
rather involved representation of unity in terms of a 
product of Dirac delta distributions on each space- 
time point, and subsequently writing these as integrals 
over auxiliary fields S a along the imaginary axis, 1 = 

j£>m n« n M s(a t s a (.x, t) - F°[s]( X , t) - o) = 

/ D[iS] J D[S ] exp [- j d d x J dt £ a S a (d, S a -F° -£»)], 
we arrive at 

(A[S])<K jD[iS] Jd[S]A[S] Jd[{] 
-J d d x J dt^S a (d,S a -F a [S]) 



exp 



exp 



-j l d x Jdt^j 



-s a c 



Performing the Gaussian integral over the noise fi- 
nally yields 



(A[S]) C = j D[S]A[S]P[S]. 
P[S]oc j D[iS ] e -* [J ' S] . 



(100) 



with the Janssen-De Dominicis response functional 



',S] = j d d x J'dt 
Yp a {d,S a - F a [S]) - S a L a [S]S a 



(101) 



The stochastic dynamics is now encoded in two dis- 
tinct sets of mesoscopic fields, namely the original 
slow variables S a and the associated auxiliary fields 
S a . Once again, the Gaussian noise normalization im- 
plies jD[iS] jD[S]e- yl[s ' s] = 1; furthermore, inte- 
grating out the auxiliary fields S a recovers the Onsager- 
Machlup functional ([99 1. 

Specifically for the purely relaxational models A and 
B, the response functional reads (we now set k%T - 1): 

W[S,S] = 

J d d x JdtJ^ [d, + D{iVf{r - V 2 )] S a 

-DS a (iV) a S a - DS a (N) a h° 
+D \ XI S a (i^yS a S l3 S l3) j , (102) 



The first two lines here represent the Gaussian action 
yio, and the term ~ u the non-linear contributions. By 
means of ( | 1 00] > and ( |102| i, the dynamical order parame- 
ter susceptibility becomes 



X aP {x-x',t-t') = 



5{S a {x,t))\ 



6hP{x', f) \h=a 
D(s a (x,t)(iV) a SP(x,t')); 



(103) 



i.e., the response function can be expressed as an expec- 
tation value that involves both order parameter S a and 
auxiliary fields S 13 , whence the latter are also referred to 
as "response" fields. Invoking causality and time inver- 
sion symmetry, it is a straightforward exercise to derive 
the fluctuation-dissipation theorem, equivalent to d85): 



X a Hx-x\t-t') = 

®( f - O ir- ( s< "( x > t)S fi tf, t')) . (104) 
of N ' 

In analogy with static, equilibrium statistical field the- 



ory (53 i, one defines the generating functional for cor- 
relation functions and cumulants, 

Z[], j] = U^Uf^f s ")\ , (105) 



4.4. Dynamic perturbation theory 

We may now proceed and treat the non-linear terms 
~ u by means of a perturbation expansion, 

n s a < s a '\ = — = — - 

ij 1=0 

Since the denominator is one owing to noise normaliza- 
tion, there are no "vacuum" contributions in this dy- 
namic field theory. Note furthermore that causality im- 
plies (S a (q, u) S 13 {q' , cij')) q - 0. From the Gaussian 
action yio (with u - 0), one immediately recovers the 



expressions ( 92 1 and ( 94 1 for Go(q, u>) and Co(q, to), re- 
spectively. Since the dynamic correlation function can 
be expressed in terms of the noise strength and the re- 
sponse function, the graphical representation in terms of 
Feynman diagrams can be based entirely on the causal 
response propagators Go(q, to), represented as directed 
lines (we use the convention that time propagates from 
right to left) that connect S 13 to S a fields, and join at ei- 
ther two-point noise or four-point non-linear relaxation 
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vertices, all subject to wave vector and frequency con- 
servation as a consequence of spatial and temporal time 
translation invariance: 



q,(0 



l 



- i(0+ D q a (r + q 2 ) 



5 «P 



2Dq"S a? 




R = _ Dq a — 
H 6 



Following standard field theory procedures, one next 
identifies the cumulants as represented by connected 
Feynman diagrams, and in complete analogy with the 
static theory establishes Dyson's equation for the full 
response propagator, G{q,toY l = Go(#, <J) _1 - £(g, <w). 
By means of the fields <X><* = 6\nZ/Sf and <$ a - 
dlnX/Sj" one constructs the generating functional for 
dynamical vertex functions via the Legendre transform 



F[(D,<5] =-]nZ[jJl 

+ J d d x J dtYjir^ + f®") , 



(107) 



iff" fT^_fy^_r[$,o]| . 

Functional derivatives then establish the following con- 
nections between two-point cumulants and vertex func- 
tions, 



r (U) ( ? ,w) = G(-q,-a))- 



^°\qM 



2Dq a 



C(q, a>) 
\G(q,co)\ 2 

Im r (U) ( ? , co) 



(108) 



(109) 



where the last relation follows from the equilibrium 
FDT (85 i. Moreover, P-°' 2 \q, at) = as a consequence 



of causality. One thus easily sees that the vertex func- 
tions are graphically represented by the one-particle 
(1PI) irreducible Feynman diagrams. 

We can now formulate the Feynman rules for the dy- 
namical perturbation expansion for the Z-th order contri- 
bution to the vertex function r (N ' N \ illustrated here for 
the one-loop graphs for F (11) and r^ 1 ' 3 -*: 

r .a 

& k v . 




1. Draw all topologically different, connected 1PI 
graphs with N out- / N incoming lines connecting 



/ relaxation vertices. Do not allow closed response 
loops (since 0(0) = in the Ito calculus). 

2. Attach wave vectors qi, frequencies w; / times ti, 
and internal indices a, to all directed lines, obeying 
"momentum- energy" conservation at each vertex. 

3. Each directed line corresponds to a response prop- 
agator Go(-q, —oj) or Goiq, ti - tj), the two-point 
vertex to the noise strength 2Dq a , and the four- 
point relaxation vertex to -Dq"u/6. Closed loops 
imply integrals over the internal wave vectors and 
frequencies or times, subject to causality con- 
straints, as well as sums over the internal vector 
indices. The residue theorem may be applied to 
evaluate frequency integrals. 

4. Multiply with -1 and the combinatorial factor 
counting all possible ways of connecting the prop- 
agators, I relaxation vertices, and k two-point ver- 
tices leading to topologically identical graphs, in- 
cluding a factor 1/7! kl originating in the expansion 
ofexp(-^ int [5,S]). 

The perturbation series then graphically becomes a 
loop expansion. For example, the propagator self- 
energy diagrams up to two-loop order are 



£ = 





With the abbreviation A{q) = Dq"{r + q 2 ), the corre- 
sponding explicit analytical expressions read: 



r (u) ( ? ,w) = ia) + Dq a 
n + 2 



r + q 



J k r + k 2 \ 6 ) J k r + 
f 1 n+2 2 f 1 

J k - (r + k' 2 ) 2 18 U J k r + k 2 



k 2 



1 



1 



X r -: k' 2 r + {q-k- k') 1 
j ito 



(110) 



io) + A(k) + A(Ar') + A(q - k - k')/ 

The renormalized noise vertex is represented by the 
vertex function F (2 ' 0) (g, a>); the first non-vanishing fluc- 
tuation correction appears at two-loop order: 
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which translates to 

T (2 '°\q,aj) = -2Dq a 



l+Dq a 



n + 2 
~18~ 



ir 



1 



1 



1 



Re 



Ik r+ k 1 J;, : r + k' 1 r + Ui - k - k') 1 
1 



ia) + A(k) + A(k') + A(q -k-k') 



(111) 



Finally, we shall require the renormalized relaxation 
vertex to one-loop order, 

T^\-3k/2;{k/2})^D(^-Ju 

i ra + 8 f _J 1 

6 Jk r + k 2 r + (q - k) 2 



1 - 



ia> + A(k) + A(q - k) 



(112) 



evaluated at equal incoming external wave vectors and 
frequencies k = (q, ai). 

4.5. Critical dynamics of the relaxational models 

We may now proceed with the perturbative renormal- 
ization of the purely relaxational models A and B, gen- 
eralizing the methods outlined in Chap. 3 to the dynami- 



cal action ( 102 1. The quadratic UV divergence (near the 



upper critical dimension d c = 4) in (110 1 is taken care 



of by an appropriate additive renormalization; since we 
are concerned with near-equilibrium kinetics here, and 
Xiq, co — 0) = x(l)> tne result is precisely the T c shift 



( 6 1 1 evaluated in the static theory. In addition to ( 63 i 
and ( |64| , we need two new multiplicative renormaliza- 
tion factors associated with the response fields and the 
relaxation rate, 



S a R = Z l J 2 S a , D R =Z D D; 



(113) 



R 



As a consequence of the FDT ( 109 », which must hold 
for both the unrenormalized and renormalized vertex 
functions, these Z factors are not independent in ther- 

/ \l/2 

mal equilibrium, but connected via Zo = (Zj /Zj) 

For model A with non-conserved order parameter 
(a = 0), extracting the UV poles in the minimal sub- 

(0, to) yields 



traction scheme from r^' 0) (0, 0) or ^F^ 



B r (U) 



n + 2 I 4 \u R 



(114) 



For model B with conserved order parameter (a = 2), 
on the other hand, the momentum dependence oc q 2 of 



the relaxation vertex implies that to all orders in the per- 
turbation expansion 

r (u) (2 = 0,co) = ico, 

d q 2^(q,oj)\ q=0 = -2D , (115) 

whence we arrive at the exact relations ZjZs = 1 and 
Zd - Zs- The conservation law thus allows us to reduce 
the dynamical multiplicative to static renormalizations. 

With ( |1 13[ ) taken into account, the renormalization 
group equation for the renormalized vertex functions 



(ji, D, t r , u R ) becomes in analogy with (|69 1 



d Ny^ + Ny s d 
^ + 2 +yDDR dD R 



d d 
+y T T R — +p u — 
ot r du R 



rf- w) = 0, (116) 



with the static RG beta function f72| and Wilson's flow 
functions ( 70 1, ( fTT) , supplemented with 



8 \ 1 Y 

y? = £f — In Zt , 
d \ D R 



7s ~ 7s 



(117) 
(118) 



owing to the FDT. Employing characteristics p —> p I to 



solve the Gell-Mann-Low RG equation ( 116 1, one has, 
in addition to ( f73j ), 



dD(€) - 
dl 



(119) 



with D( \) = D R . 

For the dynamic susceptibility near an infrared-stable 



RG fixed point, the static scaling law (75 i generalizes to 



XR (T R ,q,u>T' *p 2 i 2+y " s 



[i t D R p 2+a e 2+a+ y, 



(120) 



which allows us to identify the static critical exponents 
as before, rj = -y* s , and v = -1/y*; and in addition z = 
2 + a + y* D for the dynamic critical exponent. From the 
explicit two-loop result (114i one thus obtains for the 
0(n)-symmetric model A by inserting the Heisenberg 
fixed point ( 76 1 to order e 2 in the 4 — e expansion 



model A:z = 2 + c?7,c = 61n- - 1 + 0(e) . (121) 

Yet if the order parameter is conserved, one has y* D = 
y% , which implies the exact scaling relation 



model B : z — 4 - rj . 



(122) 
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4.6. Critical dynamics of isotropic ferromagnets 

So far we have only considered purely relaxational, 
dissipative kinetics. Often, however, the Langevin de- 
scription of critical dynamics needs to take into ac- 
count reversible systematic forces contributing to F[S] 



in (96 1. The Langevin dynamics of isotropic ferromag- 
nets provides a prominent example ||20l , (9] Chap. 6]. 
The order parameter here is a three-component vector 
field, namely the magnetization density S a (x, t), which 
represent the coarse-grained mesoscopic counterpart of 
the microscopic local Heisenberg spins, also the gener- 
ators of the rotation group 0(3). From the spin opera- 
tor commutation relation and 
Heisenberg's equation of motion, or their correspond- 
ing classical counterparts with commutators replaced by 
Poisson brackets, one readily obtains a spin precession 
term in the dynamics, in addition to the diffusive relax- 
ation of the conserved magnetization density, and con- 
served stochastic noise: 



at 



•J(x,t) : 



+DV 



8${x,f) 

(C r (x,t)^(x',t')) = -2Dk B T 
V 2 8(x - x')5(t - f')^ 



8!>(x, t) 
+ i(x,t), (123) 



(124) 



The coupled Langevin equatio ns (|123| l with the three- 
component LGW Hamiltonian ( 52 1 and noise correlator 



( 124 1 define the model J dynamic universality class. 
The associated Janssen-De Dominicis response func- 



tional comprises the model B terms, ( 102 1 with a = 2, 
and the additional reversible mode-coupling vertex, 

mj[S,S] =sj d d x jdt^ e a/3y S a 



S fi (v 2 S y +h y ) . (125) 



It is diagrammatically represented by a wave vector- 
dependent three-point vertex: 

1 



: + P 



q 

I'" 



Straightforward power counting gives [g] = jj 3 ^ 4 ^ 2 
for the mode coupling strength, which therefore be- 
comes marginal at the dynamical critical dimension 
d' c = 6, and irrelevant for d > d' c . In principle, theo- 
ries with competing upper critical dimension pose inter- 
esting non-trivial technical problems. Recall, however, 



that we are here considering near-equilibrium dynam- 
ics; thus the static critical properties completely decou- 
ple from the system's dynamics. Indeed, the dynam- 
ical critical exponent can be determined exactly from 
the underlying rotational symmetry as follows: We first 
note that according to ( 123| l an external field h y in- 
duces a rotation of the magnetization density vector: 
(s a (x, t)) h = g fidt e^y^ix, t')) h Wit). Thus, we 
obtain an exact identity linking the nonlinear suscepti- 
bility R a & = 5 2 {S a )l5hP dh^o with the linear order 
parameter response function: 



d d x ' R afiy (x,t;x-x',t -?') = 



(126) 



This identity provides crucial information for the 
renormalization of the UV singularities. In addition to 
( |1 13[ ), we define the renormalized dimensionless mode 
coupling strength 



g 2 R =Z g g 2 B dM d - 6 



(127) 



where B d = T(4 - d/2)/2 d d^l 2 . Just as for model 
B, one has to all orders in the perturbation expansion 
r (1,1) (<7 = 0, of) = iai, whence again ZjZ$ = 1. Since 



( 126 1 must hold for the renormalized susceptibilities as 



well, one then infers that Z g — Z$ . Inspection of dia- 
grams shows that the effective coupling in the dynamic 
perturbation expansion is / = g 2 /D 2 . The associated 
RG beta function becomes 



f R (d-6 + y s -2y D ) 



(128) 



Consequently, at any non-trivial, stable RG fixed point 
< /* < oo the terms in the bracket must cancel each 
other: 2 y* D — d-6 + y* s to all orders in the perturbation 
series with respect to /r. Rotation invariance and the 
conservation law thus fix the dynamic critical exponent 
in dimensions d < d' = 6 to be 



model J : z — 4 + y* D = 



d + 2-r] 



(129) 



Indeed, an explicit one-loop calculation yields y D = 
-fit + 0(u 2 R ,f^), which leads to the non-trivial model 
J RG fixed point /* = | + 0(s 2 ), where s = 6 - d. 
Note that n = for d > d c = 4, and z = 4 for 
d > d' c = 6. Since the mode-coupling vertex does 
not contribute genuinely new IR singularities, dynamic 
scaling functions for isotropic ferromagnets and related 
models can be computed to exquisit precision by means 
of the mode-coupling approximation, which essentially 
amounts to a self-consistent one-loop theory for the 
propagators, ignoring vertex corrections l20l . 
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4. 7. Driven diffusive systems 

This chapter on the application of field-theoretic 
RG methods to non-linear stochastic Langevin dynam- 
ics concludes with two paradigmatic non-equilibrium 
model systems that display generic scale invariance 
and a continuous phase transition, respectively. Both 
are driven lattice gases ED consisting of particles that 
propagate via nearest-neighbor hopping which is biased 
along a specified 'drive' direction, and is subject to an 
exclusion constraint, i.e., only at most a single parti- 
cle is allowed on each lattice site. If the system is set 
up with periodic boundary conditions, the biased diffu- 
sion generates a non-vanishing stationary mean particle 
current. At long times, the kinetics thus reaches a non- 
equilibrium steady state which turns out to be governed 
by algebraic rather than exponential temporal correla- 
tions. If in addition nearest-neighbor attractive interac- 
tions are included, the system displays a genuine non- 
equilibrium continuous phase transition in dimensions 
d > 2, from a disordered phase to an ordered state char- 
acterized by phase separation into low- and high-density 
regions, with the phase boundary oriented parallel to 
the drive and particle current. As the hopping bias van- 
ishes, the phase transition is naturally described by the 
(f-dimensional ferromagnetic equilibrium Ising model, 
since one may map the occupation numbers n, = 0, 1 to 
binary spin variables cr ; = 2« ; - 1 = +1 . 

We first consider the driven lattice gas with pure 
exclusion interactions, in one dimension also called 
"asymmetric exclusion process". In order to construct 
a coarse-grained description for the non-equilibrium 
steady state of this system of particles with conserved 
density p(x, t) and hard-core repulsion, driven along the 
'||' spatial direction on a of-dimensional lattice, one first 
writes down the continuity equation 

-S(x,i) + $ -f(x,t) = 0, (130) 
dt 

where the scalar field S(x, t) = 2p{x, t) - 1 represents a 
local magnetization in the spin language, whose mean 



remains fixed at (S (x, t)) = 0, or (p(x, t)) 



Next the 



current density must be specified; in the d ± -dimensional 
transverse sector (d ± — d - 1), one may simply assert a 
noisy diffusion current, whereas along the drive, the bias 
and exclusion are crucial: J\\ - —cD V\\S +2Dgp(l — p)+ 
(, where c measures the ratio of diffusivities parallel and 
transverse to the net current. In the comoving reference 
frame where (J\\(x, t)) = 0, therefore 



J x {x, t) = -D Vj_S(x, f) + if(x, t) , 
/||(jc,f) = -cDVyS(x,t)- -DgS(x,tf 
+«*,*), 



(131) 



with (rji) — = (£}, and the noise correlations 

(Tji(x,t)r]j(x',t')) = 

2D6(x-x')6(t-t')6 u , (132) 
<f (*, ((x, f')> = 2Dc6(x - x')5(t - t') . 

It is important to realize that Einstein's relations which 
connect the noise strengths and the relaxation rates need 
not be satisfied in the non-equilibrium steady state. 
Through straightforward rescaling of the field S(x, t), 
one may however formally enforce this connection in 
the transverse sector, say; the deviation from the Ein- 
stein relation in the parallel direction is then encoded in 
( 131 1 and ( 132 1 through the ratio < w = c/c + 1. 



The Janssen-De Dominicis response functional for 
this driven diffusive system becomes 

&[S,S] = Jd d xjdtS d,S -D(y 2 ± +cVf)S 



+Z)(vi+cVjf)S-^V||S 



(133) 



the action ( |133| > represents a "massless" field theory, 
which hence displays generic scale invariance with no 
specific tuning of any control parameters required. One 
clearly has to allow for anisotropic scaling behavior ow- 
ing to the very different dynamics parallel to the hop- 
ping bias; for example, ( |86*) l for the dynamic response 
function needs to be generalized to 

X(^M = \ q A^x{^^ d34) 

where A denotes the anisotropy exponent (A = in the 
mean-field approximation). 

Following the renormalization procedures in the pre- 
vious subsections, one first realizes that the drive gener- 
ates a three-point vertex oc g iq^\ consequently no trans- 
verse fluctuations are affected by this non-linearity, and 
= Zs = Zo = 1 to all orders in the perturbation 
expansion in g. The absence of any transverse propaga- 
tor renormalization immediately implies that 77 = and 
z = 2 in ( 134 1, which leaves merely the value of A to 
be determined. In a similar manner as for model J, one 
may in fact compute this exponent exactly; to this end, 
one observes that a generalized Galilean transformation 



S(x x ,X\\, f) — > 
iS" , X11 , t ) 



5 , (x J _,X|| 



(135) 
• Dgv t,t)-v 



leaves the Langevin equation (130i, (131 1 or equiva- 



lent action (133 1 invariant. Thus the (arbitrary) speed 



v must scale as the order parameter field S , and since 
Z D = 1 = Z s , neither can the coupling g be altered by 
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fluctuations, or ( |135| l would be violated for the renor- 
malized theory. Hence Z g — 1 as well, and the only 
remaining non-trivial renormalizations are those for the 
dimensionless parameters c R -Z c c and c R = Z- c c. 

An explicit one-loop calculation establishes the exis- 
tence of an IR-stable RG fixed point for the coupling 



v = g 2 /c 3/2 , v R =Z^vC di X 



3/2, 



4-2 



(136) 



only be retained in the transverse sector. For this driven 
model B, the effective critical action thus becomes 

M[S,S]= Jd d xjdts\d,S - DV\ (r-V 2 x )s 

-DcV\S +d{v 2 J - | VyS 2 - ~ V^S 3 ) , 

(141) 



with C d = F(2 - d/2)/2 d - l n d/2 , and identifying d c = 2 
as the upper critical dimension for this problem. Eval- 
uating the one-loop fluctuation corrections to the longi- 
tudinal propagator, one finds 



and (134 1 is further generalized by adding a relevant 
temperature variable (now A = 1 in mean-field theory): 



X(T±,q±,q\\,u) = \q±\ 
T 1\\ 



-2+1 



7c = -T7 ( 3 + W R> » 

to 



(137) 



X 



\qA llv \q±\ 1+K \qJ 



(142) 



for the anomalous scaling dimensions of c and c, or 



J3 W = w R (y e - y c ) 

= ~^§( w * ~ - 3) 



B v = v R \d-2- -y t 



(138) 
(139) 



Straightforward power counting for the non-linear 
couplings yields [g 2 ] = fi 5 ~ d and [u] = therefore 
the upper critical dimension is raised to d c = 5 (com- 
pared to both the non-critical driven lattice gas and the 
equilibrium model B), and fluctuation corrections are 
dominated by the drive, while the static coupling u is 
(dangerously) irrelevant near d c . The three-point vertex 
oc g iq\\ again does not allow any renormalizations in the 
transverse sector, whence Zj = Z $ = Z p = 1; conse- 



quently 77 = 0, v = L and z — 4 in ( 142 1 to all orders in 



for the associated RG beta functions of the ratio w — c/c 
and the non-linear coupling v. At any non-trivial RG 



fixed point < v* < 00, ( 138 1 implies that either w* N = 3 



the perturbation series, leaving only the anisotropy ex- 
ponent to be determined. As before Z g = 1 follows from 
Galilean invariance, imposing a simple structure for the 



or w* E = 1, but the latter is obviously stable; in the RG beta function for the effective coupling ( 136 1: 
asymptotic scale-invariant regime, the Einstein relation 
is evidently restored, and the system effectively equili- 
brated. Moreover, in dimensions d < d c = 2, ( |139| > 
leads to the exact scaling exponents 



p v = v R \d-5- - y t 



(143) 



2^ = 2 ~ d 
' 2 3 



zil 



1 + A 5-d 



(140) 



In dimensions d < d c — 5, at any non-trivial and finite 
RG fixed point, the scaling exponents are thus forced to 
assume the values 



At d = 1, specifically, one has z\\ - |, which cap- 
tures the dynamic scaling for the asymmetric exclusion 
process. In one dimension, the driven lattice gas with 
exclusion in fact maps onto the noisy Burgers equa- 
tion for equilibrium fluid hydrodynamics, and also to 
the Kardar-Parisi-Zhang equation for curvature -driven 
surface or interface growth. 

We finally briefly summarize the RG analysis for the 
driven lattice gas with conserved total density and at- 
tractive Ising interactions between the particles (and 
"holes") at its critical point. Since the system orders 
in stripes along the drive direction, only the transverse 
fluctuations become critical. Therefore one must amend 



1 



S-d 



12 



i + A n-d 



• (144) 



the response functional ( 133 1 with a higher-order gra- 
dient term and non-linearity akin to the scalar model 



B, see ( 102 1 with a — 2; yet the noise terms too need 



These last examples clearly demonstrate how the pow- 
erful field-theoretic RG approach can help to exploit the 
basic symmetries for a given problem, allowing to de- 
termine certain non-trivial scaling exponents exactly. 

5. Scale Invariance in Interacting Particle Systems 

This last chapter details how the stochastic kinetics of 
classical interacting (reacting) particle systems, defined 
through a microscopic master equation, can also be 
mapped onto a dynamical field theory in the continuum 
limit |22|-[24|. For at most binary reactions, one can 
thus derive a corresponding mescoscopic Langevin rep- 
resentation, typically with multiplicative noise terms. 
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Furthermore, RG tools may be applied to extract the 
infrared properties in scale-invariant systems, as will 
be exemplified for diffusion-limited annihilation pro- 
cesses ll22ll24l[T9l , and for non-equilibrium phase tran- 
sitions from active to inactive, absorbing states, where 
all stochastic fluctuations cease |24| |27|. Stochastic 
models in population dynamics and ecology are natu- 
rally formulated in a chemical reaction language, and 
hence amenable to these field-theoretic tools ||28l [291 . 

5.1. Chemical reactions and population dynamics 

Let us thus consider (classical) particles of various 
species A, B, ... on a cZ-dimensional lattice that propa- 
gate by hops to nearest-neighbor sites, and either spon- 
taneously decay or produce offspring, and/or upon en- 
counter with other particles, undergo certain "chemical" 
reactions with prescribed rates. Our goal is to systemati- 
cally construct a continuum description of such stochas- 
tic particle systems that however faithfully encodes the 
associated intrinsic reaction noise, and consequently al- 
lows us to properly address the effects of statistical fluc- 
tuations and spatio-temporal correlations [22, 24 HT9ll . 

To set the stage, we introduce three characteristic ex- 
amples that we shall explore in more detail below. First, 
we address the general single-species irreversible anni- 
hilation reaction kA — > mA with integers m < k, and 
rate A k . Assuming the system to be well-mixed, one 
may neglect spatial variations and focus on the mean 
particle density a(t) = (a(x,t)). Ignoring in addition 
any non-trivial correlations, one can write down the rate 
equation for this stochastic process, which in essence 
thus constitutes the simplest mean-field approximation: 



consequence, one obtains slower decay power laws than 



daif) 
~dt~ 



- -(k - m) A k a(tf 



(145) 



For k — 1 (and m = 0), this just describes spon taneo us 
exponential decay, a(t) = a(0)e~ A ' '; for k > 2, ( 145 1 is 
easily integrated with the result 

a(t) = [a(0) l - k + (k- m)(k - 1) A k f ]" 1/(fc_1) ', (146) 

For the kth order annihilation reaction, the particle den- 
sity decays algebraically a(t) ~ (Ak f) _1 ^* _1) at long 
times t » zlr 1 , with an amplitude that does not even de- 
pend on the initial density a(0) anymore. The replace- 
ment of an exponential decay by a power law signals 
scale invariance and indicates the potential importance 
of fluctuations and correlations. Indeed, the annihilation 
kinetics generates particle anti-correlations, whence the 
long-time kinetics is dominated by the ensuing deple- 
tion zones in low dimensions d < d c (k) that need to 
be traversed by any potentially reacting particles. As a 



predicted by the mean-field rate equation ( 146 1. 

Next, we allow competing reactions, namely decay 
A — > (the empty state) with rate k, and the reversible 
process A ^ A+A with forward / backward rates cr and 
A, respectively. Again, we begin with an analysis of the 
rate equation for this set of reactions, 



= (cr- K )a{t)- Aa(t) 2 

at 



(147) 



which obviously predicts a continuous non- equilibrium 
phase transition at cr c — k: For cr > k, the mean par- 
ticle density approaches a finite value, a(t — > oo) — > 
aca - {cr - k)/A. One refers to this state as an ac- 
tive phase; ongoing reactions cause the particle number 
to fluctuate about its average. On the other hand, for 
cr < k, the density can only decrease, whence ultimately 
a(t) — > 0; in this inactive phase, all reaction processes 
terminate since they all require the presence of a par- 
ticle. Such a state is therefore called absorbing: once 
reached, the stochastic dynamics cannot escape from it 
anymore. Right at the critical point cr = k, one recov- 
ers the long-time algebraic decay of the pair annihilation 
process, a(t) ~ (A f) , this suggests the interpretation of 



( 146 1 as a critical power law induced by the precise can- 



cellation of the contributions from first-order reactions 
that enter linearly proportional to the particle concen- 
tration. The obvious issues to be addressed by a more 
refined theoretical treatment are: How can internal re- 
action noise and correlations be systematically incorpo- 
rated? What is the upper critical dimension d c below 
which fluctuations crucially alter the mean-field power 
laws? Can certain universality classes be identified, and 
the associated critical exponents be computed, at least 
perturbatively in a dimensional expansion near d c l 

Finally, let us address a prominent textbook ex- 
ample from population dynamics, namely the classic 
Lotka-Volterra predator-prey competition. Invoking the 
stochastic chemical reaction framework, this model is 
defined via spontaneous death A — > (rate k) and birth 
B — > B + B (rate cr) processes for the "predators" A 
and "prey" B; absent any interactions between these two 
species, the predators must go extinct, while the prey 
population explodes exponentially. Interesting species 
competition and potentially coexistence is created by 
the binary predation reaction A + B — > A+A (with 
rate A). The associated coupled rate equations for the 
presumed homogeneous population densities read 

8a(t) 



dt 
db(t) 



Aa(t)b(t) - Ka(t) 



= crb(t) - Aa(t)b(t) . 



(148) 
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In this mean-field approximation, one easily confirms 
the existence of a conserved first integral for the ordi- 
nary differential equations ( |148[ ): The quantity K(t) = 
A[a(t) + b(t)] - ln[a(t) cr b(t) 1 '] = K(0) remains unchanged 
under the temporal evolution. Consequently, the mean- 
field trajectories are closed orbits in the phase space 
spanned by the population densities, and the dynamics 
is characterized by regular population oscillations, de- 
termined by the initial state. This is clearly not a bio- 
logically realistic feature, and indeed represents an arti- 
fact of the implicit mean-field factorization for the non- 
linear predation process. Upon including the internal 
reaction noise and spatial degrees of freedom with dif- 
fusively spreading particles, as, e.g., in individual-based 
Monte Carlo simulations, one in fact observes striking 
"pursuit and evasion " waves in the species coexistence 
phase that generate complex dynamical patterns, locally 
discernible as erratic population oscillations which ul- 
timately become overdamped in finite systems. More- 
over, if the local "carrying capacity" is finite, i.e., only a 
certain maximum number of particles may occupy each 
lattice site, there emerges a predator extinction thresh- 
old which indicates a continuous phase transition to 
an absorbing state, namely the lattice filled with prey. 
Stochastic fluctuations as well as reaction-induced noise 
and correlations are thus crucial ingredients to properly 
describe the large-scale features of spatially extended 
Lotka-Volterra systems even and especially far away 
from the extinction threshold (for recent overviews, see 

Refs. nana). 

5.2. Coherent- state path integral for master equations 
In the following, reacting particle systems on a d- 
dimensional lattice shall be defined through the associ- 
ated chemical master equation governing a Markovian 
stochastic process with prescribed, time-independent 
transition rates. Any possible configuration at time t 
of the stochastic dynamics is then uniquely character- 
ized by a list of the integer occupation numbers n, 
, 1,2.... for each particle species at sites i. The mas- 
ter equation governs the temporal evolution of the cor- 
responding probability distribution P({«,}; t) through a 
balance of gain and loss terms induced by the reac- 
tion processes; for example, for the binary reactions 
A + A -> and A + A -» A with rates A and A': 
d 

P(nu t) - A (ra ; + 2) (n,- + 1) P(. . . ,m + 2, . . . ; t) 



dt 



+A' (rii + 1) m P(. . . , Hi + 1, . . . ; t) (149) 
-(A + A')n i (n i -l)P(...,n h ...;t), 



with, say, an uncorrected initial Poisson distribution 
P({n i },0) = Ui{K le ~' k /"4 



Since the dynamics merely consists of increasing or 
decreasing the particle occupation numbers on each site, 
it calls for a representation through second-quantized 
bosonic ladder operators, at least if arbitrary many par- 
ticles are permitted per site, with standard commuta- 
tion relations [a,-,o/] = 0, [a,-, at] = <Sy, and an empty 
vacuum state |0> that is annihilated by all operators a,, 
a,|0) = 0. The Fock space of states |{«,}) with n,- 
particles on sites i is then constructed through multi- 
ple creation operators acting on the vacuum, |{«,}) = 
n,(flj) |0); note that a different normalization from 
standard many-particle quantum mechanics has been 
implemented here. Thus, a-, |n,) = n, |n, - 1) and 
a\ \rii) = \rti + 1), whence the states |{«,}) are eigen- 
states of «,■ = a.'fl/ with eigenvalues «,-. Next one defines 
the formal state vector |<I>(f)) = HinA P({ n iY> W n i))i 
whose temporal evolution is determined by the mas- 



ter equation (149 1, and may be written in terms of a 
time-independent quasi-Hamiltonian or Liouvillian H 
that can be decomposed into a sum of local operators: 



dt 



\<W> = -H |<D(f)> ,H=y H t {a], ad . (150) 



Note that ( 150 1 constitutes a non-Hermitian imaginary- 
time Schrddinger equation, with the formal solution 
|O(0> = exp(-^)|3>(0)>. 

For example, the quasi-Hamiltonian in this Doi- 
Peliti bosonic operator formulation reads for diffusion- 
limited annihilation and coagulation reactions 

ff = £ 2(4 -«})("<•-";) ( 151 ) 

«7> 

-EK 1 - fl 'V +A '( 1 - fl i) a X . 

where the first line represents nearest-neighbor hopping, 
and the second encodes the processes in ( 149 1. For 



each stochastic reaction, H contains two contributions: 
the first one directly reflects the physical process un- 
der considerations, i.e., annihilation and production of 
particles, whereas the second term carries information 
on the order of the reaction (which powers of the par- 
ticle concentrations enter the rate equations). In order 
to access the desired statistical averages with the time- 
dependent probability distribution f({n,};f), one needs 
the projection state (P\ = <0| Hi e '", with <P|0> = 1; the 
mean value for any observable F, necessarily a function 
of all occupation numbers n,-, at time t then follows as 



<F(0> = E F({/Ii})p({rei};f) 

{",) 

= {'P\F{{a]a i })m))- 



(152) 
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Probability conservation implies that 1 = ( , P|<E>(?)) = 
(f>| e - fff |O(0)>, and therefore (P\H = 0. By means 
of Je a ,a'] = e", one may commute the product e^ ,a < 
through the quasi-Hamiltonian H, which effectively re- 
sults in the operator shifts a! — > 1 + a! . H must conse- 
quently vanish if all creation operators are replaced with 
l,Ht(a. — > 1 , a,) = 0. In averages, one may thus also re- 
placea.'a, — > a,; e.g., for the particle density one obtains 
simply a(t) = (a,), while the two-point occupation num- 
ber operator product becomes a\ a,- qlfl j — > a-fiij + a,-ay. 
Starting with the Hamiltonian ( 150 1, and based on the 



expectation values ( 152 1, one may invoke standard pro- 



cedures from quantum many -particle theory to construct 
a path integral representation based on coherent states, 
defined as eigenstates of the annihilation operators a, 
with arbitrary complex eigenvalues ip,: a t = (p, \<pj). 
It is straightforward to confirm 

|^> = e X p(-i|^| 2 + 0i4)|O>, (153) 



\mxm\. 



(154) 



The closure relation ( |154[ ) demonstrates that the coher- 
ent states for each site i form an overcomplete basis of 
Fock space. Upon splitting the time evolution into in- 



finitesimal steps, and inserting ( 154 1 into ( 152 1 with the 
formal solution of (|150[), one eventually arrives at 



(F(t)} kJYI D[< ^ i] F({(/>i}) e~ m *\ 

i 

f dt[4>'t 3 ( 0, + //(0;,0,)]-«o0*(O) 



f 

Jo 



In the end we take the continuum limit 0,(f) — > a Q 
»/(x,t) (with lattice constant a^), and </>*(f) — > {jj{x,t); 
for diffusively propagating particles, the ensuing "bulk" 
action becomes 



iff] = J d d x J f dt (d, - D V 2 ) (A 



(156) 



Here, 1ri r denotes the contributions stemming from the 
stochastic reaction kinetics; e.g., for pair annihilation 
and coagulation, 

= -4} - - ^'(i - <a)<a<a 2 • (157) 

Appropriate factors of ao were absorbed into the con- 
tinuum diffusion constant D and reaction rates A,A' . It 



is worthwhile emphasizing that the actions ( 155 i based 
on a master equation should be viewed as microscopic 
stochastic field theories, which may well require addi- 
tional coarse-graining steps. Yet the internal stochastic 
dynamics of the master equation is faithfully and con- 
sistently accounted for, since aside from the continuum 
limit no approximations have been invoked; specifically, 
no assumptions on the form or strength of any noise 
terms have been made. As exemplified next for the 



action ( 157 1, the Doi-Peliti coherent-state path integral 
representation of stochastic master equations may serve 
as a convenient starting point for systematic analytical 
approaches such as field-theoretic RG studies. 

5.3. Diffusion-limited annihilation processes 

Pair annihilation A +A — > or coagulation A +A — > A 
represent the perhaps simplest but non-trivial diffusion- 
limited reactions. In order to reach beyond the mean- 



field rate equation approximation ( 146 1, we explore the 



corresponding Doi-Peliti field theory ( 156 1 with the 



specific reaction Hamiltonian ( 157 1. First we note that 
the associated classical field equations — = 

SJl/dij/ are solved by {(/ — 1 (which just reflects proba- 
bility conservation) and 



difi(x, t) 
dt 



D V iff(x, t) - (2 A + A') ifr(x, t) , (158) 



i.e., the rate equation for the local density field xpix, t) 
augmented by diffusive spreading. It is convenient to 
shift the conjugate field about the mean-field solution, 
\j/{x, f) = l+ iff(x, t), which turns the reactive action into 

'Hrty, iff) = (2A + A')^if, 2 + (A + A') i// 2 i/r 2 . (159) 

Since the annihilation and coagulation processes oc A, A' 
generate the very same vertices, we conclude that aside 
from non-universal amplitudes, both diffusion-limited 
reactions should follow identical scaling behavior. One 
may also formally interpret the ensuing field theory as a 
Janssen-De Dominicis response functional (JTOTJ orig- 
inating from a "Langevin equation" ( 158 i with added 
white noise £(x, f), whose second moment ( |97] > is given 
by the functional L[if/] = -(A + A') xjj 2 < 0. This negative 
variance, which reflects the emerging anf/-correlations 
for surviving particles that are induced by the annihila- 
tion reactions, of course implies that a Langevin repre- 
sentation is not truly feasible for this stochastic process. 
One must also keep in mind that the fields ijj and are 
complex-valued; indeed, the reaction noise can be recast 
as "imaginary " multiplicative noise oc iif/(x, i) £(x, t) in 
the associated stochastic differential equation. 
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As coagulation thus falls into the same universality 
class as annihilation, let us more generally study the 
Doi-Peliti action for k-particle annihilation k A — > 0, 



(160) 



The corresponding mean-field rate equation ( |145| l pre- 
dicts algebraic decay a(t) ~ (A k t)~ l ^ k ~ l) at long times. 
Since the field \fi appears to higher than quadratic power 
for k > 3, no (obvious) equivalent Langevin descrip- 
tion is possible for triplet and higher-order annihilation 
reactions. With [\j/(x, t)] = 1 and [i/f(x, f)] = H , as 
(ifr(x, f)) = a(t) is just the particle density, power count- 
ing gives [A k ] = ii 2 ~( k ~ l ^ d ; the upper critical dimension 
thus is d c (k) = 2/(k - 1) for kth order annihilation, and 



one expects the mean-field power laws ( 146 1 to be ac- 
curate for k > 3 in all physical dimensions d > 1; for 
triplet reactions, one should encounter merely logarith- 
mic corrections in one dimension, where genuine non- 
trivial exponents ensue only for pair annihilation. 
Even for k — 2, one cannot construct any loop graphs 



from the vertices in ( 160 1 that would modify the mass- 
less diffusion propagator, implying that 77 = and z — 2. 
This leaves the task to determine the reaction vertex 
renormalization, which can also be achieved to all or- 
ders by summing the diagrammatic geometric series 
(essentially a Bethe-Salpeter equation; here for k = 3): 

+ 



+ 



With the factor B kd = k\ T(2- d/d c ) d c /k d/2 (An) dld % one 
finds the renormalized reaction rate 



A 

8r - Z g — B M p 



Z g ' = 1 + 



2(l-d/d c ) 

AB kdi r^- d ^ 



(161) 



D (d c - d) 

and the exact RG beta function and stable fixed point 

2gR 



(d - d c + g R ) , g 



(162) 



Next we write down the Gell-Mann-Low RG equa- 
tion for the particle density a(t), applying the matching 
condition (jjt) 2 = 1 /Dt : 



d(D t) dn dg R 

a(jx,D,n ,g R ,f) = , 



with the solution 

a(jx,D,n ,gR,t) = 

(D M 2 ty d > 2 a(n (Dp 2 t) d/2 ,g(t)) . (164) 

The particle density at time t naturally depends on its 
initial value «o, clearly a relevant parameter in the RG 
sense. One therefore needs to establish through explicit 
calculation that the (tree level) scaling function a re- 
mains finite to all orders as «o — * 00 ■ In the end, ( |164| i 
yields for pair annihilation, 

k-2: d<2: a(t) ~ (D iy dl2 , 

d-2: a(t) ~ (D ln(D f) , (165) 

d > 2 : a(t) ~ (At)' 1 ; 

while for the triplet reaction 

k = 3 : d = 1 : a(t) ~ [(Dt)' 1 ln(Df)] 1/2 , 

d> 1 : a(t) ~ (Aty l/2 . (166) 

At low dimensions d < d c (k) = 2/(k — 1), the den- 
sity decay is slowed down as compared to the mean- 
field power laws by the emergence of depletion zones 
around the surviving particles. Further annihilations 
require that the reactants traverse the diffusion length 
L(t) ~ (Dt) 1 ^ 2 , which sets the typical separation scale; 
the corresponding density must then scale as L(t)~ d , Be- 
yond the upper critical dimension d c (k), the system be- 
comes effectively well-mixed, diffusion plays no limit- 
ing role, and the time scale is set by the reaction rate. 

5.4. Phase transitions from active to absorbing states 

Turning to our second example in the introductory re- 
marks, we now investigate diffusing particles subject to 
the competing reactions A — > and A ^ A + A; adding 
the diffusion term to ( |147| i, we arrive at the rate equation 
for the local particle density, 

da(x, t) 



dt 



= —D{r— V 2 ) a(x, f ) - A a(x, t) 2 



(167) 



where r = (k-ct)/D; in mathematical biology and ecol- 



(163) 



ogy, the partial differential equation ( 167 1 is known as 
the Fisher-Kolmogorov equation, and for example has 
been used to study population invasion fronts into empty 
regions. We shall instead focus on the critical region 
where the control parameter r — > 0, and a continuous 
non- equilibrium phase transition from an active to an 
inactive and absorbing state occurs. 

The Doi-Peliti field theory action ( 156| l capturing the 
above reactions reads 

$0,ip] = j ddx J dt ${d,-DV 2 )<p (168) 
-k(i - tff)i// + oil - iA)^- A(l - <a)*A *A 2 
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Upon shifting and rescaling the fields according to 
ft(x,f) = 1 + y/A/o- S (x, t) and ij/{x, t) = y/cr/AS(x, t), 
one arrives at 

W[S,S] = jd d xjdt(s[d,+D(r-V 2 )]s 

-u{s -S)SS +AS 2 S 2) j, (169) 

where the three-point vertices have been symmetrized 
and now are proportional to the coupling u = Vcr A with 
scaling dimension [u] = p?^ d ^ 2 . The associated upper 
critical dimension is therefore d c = 4, and the annihila- 
tion four-point vertex [A] = p 2 ~ d consequently is irrel- 
evant in the RG sense near d c . In the effective critical 
action, one may set A — > 0, whereupon ( |169| l reduces to 
the familiar Reggeon field theory action, which is invari- 
ant under rapidity inversion S (x, t) <-> -S (x, -t). View- 



ing ( 169 1 as a Janssen-De Dominicis functional (101 



it is equivalent to the Langevin equation that amends the 



Fisher-Kolmogorov equation ( 167 i with a noise term, 
dS (x, t) 



dt 



= D(y 2 -r)S(x, i)~uS(x, t) 2 +((x, t) , (170) 



with (£(x, t)) = and the multiplicative "square-root" 
noise correlator 

tX(x' , 0> = 2u S (x, f) S(x - x')6(t - O . (171) 

Drawing a space-time plot (time running from right 
to left) for the branching A -> A + A, death A — > 0, 
and coagulation A + A — > A processes, starting from 
a single occupied site, as depicted below, one realizes 
that they generate a directed percolation (DP) cluster, 
with "time" playing the role of a specified "growth" di- 
rection. The field theory (169 1 (with A — 0) thus also 
describes the universal scaling properties of critical DP. 




Indeed, one expects that active to absorbing state 
phase transitions should generically be captured by this 
DP universality class, namely in the absence of coupling 
to other slow conserved modes, disorder, or special ad- 
ditional symmetries. The origin for this remarkable DP 
conjecture becomes evident in a complementary coarse- 
grained phenomenological approach that will be framed 
in the language of epidemic spreading |27|. Consider 
the following simple epidemic process: 



1. A "susceptible" medium is locally "infected", de- 
pending on the density of "sick" neighbors. In- 
fected regions may later recover. 

2. The "disease" extinction state is absorbing. 

3. The disease spreads diffusively via infection, see 1. 

4. Other fast microscopic degrees of freedom are in- 
corporated as random noise. Yet according to con- 
dition 2, noise alone cannot regenerate the disease. 

These decisive features can be encoded in a meso- 
scopic Langevin stochastic differential equation for the 
local density n(x, t) of "active" (infected) individuals, 

dn(x, t) 



dt 



D (V 2 - R[n(x, t)]) n(x, t)+£(x, t) , (172) 



with the reactive functional R[n], {((x, t)) = 0, and noise 
correlator L[n] = n N[n]. In the spirit of Landau theory, 
near extinction one may expand these functionals in a 
Taylor series for small densities, 



r x : R[n] = r+u «+..., N[n] = v+. 



(173) 



where higher-order terms are irrelevant in the RG sense. 
After rescaling, the corresponding Janssen-De Domini- 



cis response functional ( 101 1 becomes identical to the 
Reggeon field theory action. 

We now proceed to analyze the dynamic perturba- 



tion theory and renormalization for the DP action ( 169 1 
to one-loop order. The only singular vertex functions 
are the propagator self-energy r (1,1) (g, ai) and the three- 
point functions r (12) = -F (2,1) , owing to rapidity inver- 
sion symmetry, with the lowest-order Feynman graphs: 





Explicit evaluation of the self-energy yields 

r (u, (<7, w) = iu> + D(r + q 2 ) 
.2 r ! 



u z r 



iu/2D + r + q 2 /4 + k 2 



(174) 



As in critical statics, one needs to first ensure the criti- 
cality condition, namely that r (1,1) (0, 0) = at the true 



percolation threshold r = r c . To one-loop order, ( 174 1 
results in the shift (additive renormalization) 



+ 0(u 4 ) . 



(175) 



f 1 

rc ~ DiJ k r c+ k 2 

and inserting t = r - r c in ( 174| > subsequently leads to 

r (u) ( ? ,w) = iu) + D(r + q 2 ) (176) 

u 2 C iu/2D + T + q 2 /4 
~Dj k 



k 2 (Uo/ID + t + q 2 /4 + k 2 ) 
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The three-point vertex function at vanishing external 
wave vectors and frequencies finally becomes 

2m 2 



r (U) ({0},({0}) = -2 M 



1 



f— 

Jk (t + k 



.(177) 



D 2 J k ( T + k 2 ) z 
For the multiplicative renormalizations, we follow 



(63 1 and d64l, but with the sligh t modification u R 



Z u uAJ pF~^' , l 2 , as well as ( 1 1 3 » with Zj = Zs ow 



ing to rapidity inversion invariance. Thus one obtains 
Wilson's RG flow functions to one-loop order, 

VR 



ys 



j + 0(v 2 R ) 



JD 



+ OK) . 



Jr = -2 + ^ + 0(V 2 R ) , 



(178) 



with the effective coupling and associated beta function 

u 2 



Z 2 - 2 



v R =-,- 2 A df < 



jl-4 



-e + 3v R + 0(v R ) 



(179) 



(180) 



Below the upper critical dimension d c = 4, the IR-stable 
RG fixed point 



"dp = 3 + 0(e z ) 



(181) 



appears, and solving the RG equation for the two-point 
correlation function in its vicinity yields 

1 ~ q 2 & 

q oj \ 



C R (T R ,q,u>) 

C R (r R V .V . 



(182) 



' pt D R p 2 ( 2 +v"D> 

This allows us to identify the three independent critical 
exponents for directed percolation to order e = 4 - d: 

n = -7s =-^+0(e 2 ), 



2 + 7 D 



+ 0{e 2 ) 



(183) 



e 

12 



0(e 2 ) 



The DP universality class also applies to many ac- 
tive to absorbing state phase transitions with more than 
just one particle species. As an example, consider the 
predator extinction threshold in the spatially extended 
stochastic two-species Lotka-Volterra model with finite 
carrying capacity discussed in the chapter introduction. 
The associated Doi-Peliti field theory action reads 

S [a, a; b, b] = J~d d x J"dt a [d t — D^V 2 ) a 

+K(a-\)a + b(d,-D B V 2 )b (184) 

+<x(l -b)bbe- p '' hb + A[b-a)aab , 



where diffusive spreading has been assumed, and the ex- 
ponential term in the prey production term takes into 
account the local restriction to a maximum particle den- 
sity p. As usual, one applies the field shifts a = 1+5, 
b — 1 + b; realizing that p~' constitutes an irrelevant 
perturbation (since the density scales as [p] = p d ), we 
furthermore expand to lowest order inp -1 , which yields 

S [5, a; b, b] = Jd d x Jdt a (d, -D A V 2 + Kja 
+b{d,-D B V 2 -o)b-o-b 2 b (185) 
+o-p"'(l +bf~bb 2 -A.{\ +a)(a-b)ab 



Near the predator extinction threshold, the prey almost 
fill the entire system. We therefore define the properly 
fluctuating fields c — b s — b with b s x p and (c) = 0, 
and c = —b. Rescaling to (p — yfcrc and (f> = -\fcrc, 
and noting that asymptotically <x — > oo under the RG 
flow since [cr] = p 2 , the ensuing action simplifies dras- 
tically. At last, we add the growth-limiting reaction 
A + A — > A (with rate r); the fields (p and (f> can then 



be integrated out, leaving Reggeon field theory ( 169 1 as 
the resulting effective action, with the non-linear cou- 
pling u = ■\frTb~ s . 



6. Concluding Remarks 

These lecture notes can of course only provide a very 
sketchy and vastly incomplete introduction to the use of 
field theory tools and applications of the renormaliza- 
tion group in statistical physics. I have merely focused 
on continuous phase transitions in equilibrium, dynamic 
critical phenomena in simple relaxational models, and 
a few examples for universal scaling behavior in non- 
equilibrium dynamical systems. Among the many top- 
ics not covered or even mentioned here are critical phe- 
nomena in finite and disordered systems; universality 
classes of critical dynamics with reversible couplings to 
other conserved modes; universal short-time and non- 
equilibrium relaxation scaling properties in the aging 
regime; depinning transitions and driven interfaces in 
disordered media; spin glasses and structural glasses; 
and of course phase transitions and generic scale in- 
variance in quantum systems. Nor have I addressed 
powerful representations through supersymmetric or 
conformally invariant quantum field theories, Monte 
Carlo algorithms, or numerical non-perturbative RG ap- 
proaches, since the latter will be covered elsewhere in 
this volume; for their applications to non-equilibrium 
systems, see, e.g., Ref. QUI . 
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Nevertheless, I hope to have conveyed the message 
that methods from field theory are ubiquitous in statis- 
tical physics, and the renormalization group has served 
as a remarkably powerful mathematical tool to address 
at least the universal scaling aspects of cooperative be- 
havior governed by strong correlations and fluctuations. 
Thus, the RG has become a cornerstone of our under- 
standing of complex interacting many-particle systems, 
and its language and basic philosophy now pervade the 
entire field, with applications that increasingly reach out 
beyond fundamental physics to material science, chem- 
istry, biology, ecology, and even sociology. 

Finally, I would like to express my gratitude to the 
organizers (and their funding agencies) for the opportu- 
nity to attend and lecture at the 49th Schladming The- 
oretical Physics Winter School. I thoroughly enjoyed 
the stimulating and informal atmosphere in the Styr- 
ian Alps, and profited from many discussions with col- 
leagues and students. I can only hope that all attendants 
learned as much about new and exciting developments 
in theoretical physics from my fellow lecturers as I did. 
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